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ABSTRACT 


Three lower confidence interval estimation procedures for system reliability of co- 
herent systems with cyclic components are developed and their accurecy evaluated by 
Monte Carlo methods. Each method uses estimates of the ratios of component unreli- 
abilities and the Poisson approximation to the binomial distribution to obtain the 
equation for the lower confidence limit. This is an extension of a method previously 
reported in the literature which has been.shown to be fairly robust. The procedures de- ae 
*veloped here can be combined with similar procedures already developed for systems . 

with continuous components. The combined procedure may yield a reasonably accurate 
lower confidence interval procedure for the reliability of coherent systems with mixtures 


of continuous and cyclic components. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs devcioned in this research may not 
have been exercised for all cases of interest. While every effort has becn made, within 


the Ume available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. BACKGROUND 


A coherent system is any system whose reliability is not reduced when the reliability 
of any of its components is increased. This thesis deals with coherent systems that have 
cyclic components. Cyclic components perform one or more repetitions of a single 
function during their mission. Examples are components that have power-turn-ons and 
power-turn-offs, switches and hydraulic units. Any component that performs a function 
which is regarded as a success or failure is a cyclic component. Consequently a contin- 
uously operating electronic component becomes a cyclic component if the only data for 
that component is the number of successes and failures in repeated trials. The number 
of trials to first failure of cyclic components have discrete probability distributions. 
Numerous approximate interval estimation procedures have been developed for system 
reliability of coherent systems with cyclic components. Such a method developed by 
Myhre and Saunders [Ref. 1: p.37] uses a likelihood ratio method. C. R. Rao [Ref. 2 } 
developed a maximum likelihood method and Easterling [Ref. 3 ] proposed a modified 
maximum likelihood method. For special structures, such as series systems, exact 
methods were developed by Winterbottom [Ref. 4: pp-782-787]. Other approximate and “™" 
asymptotic methods were developed by Mann and Grubbs [Ref. 5 ]. The accuracy of 
some of these approximate procedures has been studied for specific sample size cases for 
structures of order two or three by Mann and others [Ref. 6] and Winterbottom [Ref. 
4: pp.782-787]. Woods and Borsting [Ref. 7 ] developed an approximate procedure 
which was modified extensively by Mann and Grubbs [Ref. 5: pp.335-347]. None of 
these discrete interval estimation procedures cited above can be readily used in con- 
junction with continuous component data to obtain interval estimates for the reliability 
of complex systems that have mixtures of cyclic and continuous components. 

In this thesis we attempt to establish feasible interval estimation methods for the 
reliability of coherent systems with cyclic components. These methods have a common 
feature that allow them to be combined with similar methods that use continuous data. 
The combination of these methods may provide interval estimates for the reliability of 
systems with cyclic and continuous components. This method is an extension of a 
method developed by Lomnicki [Ref. 8 ] and extended by Myhre and others [Ref. 9: 
p.213}. They compare the accuracy of their procedure with the accuracy of three of the 


procedures cited above and provide bounds for error between the true lower confidence 














bound and their approximate lower bound. These bounds are used to show that their 
procedure is fairly robust against errors in assumptions about the ratios between com- 
ponent unreliabilitics which are necded to use their procedure. The procedure developed 
in this thesis attempts to capitalize on this robustness to extend their method by esti- 
mating these ratios from the test data . Specifically, let R, denote the reliability of com- 
ponent i and denote component unreliability by ¢g.=1—R, Suppose the component 
unreliabilities‘of a coherent system of size k are denoted by q,, 9,,...,9,. Let 


Im = kale Jas 9 Wp} (1.1) 


qi Squeys ; : 
and a,= a Then the reliability function h(q,, 9, ... ,9,), of the system can be written 
aS N(qny @, Gy, ...,,)- Since the system is coherent, a lower confidence limit on system 
reliability Re is given by 


mat: ore Rsp = h(@.,. Ayy very ay) (J sD} 


where q, is the corresponding upper confidence limit on q,. In this thesis the a, are es- 
timated from the data. The method developed by Myhre and others [Ref. 9: pp.216-223] 
assume the a, are known. 

Mvhre and others [Ref. 9: p.223] also apply this failure rate ratio concept to systems 
whose components have exponentially distributed failure times with failure rates 4, In 


this case, 
Am = Max{A,, 42, -- Ax} (1.3) 
se 1.4 
a= hn ( . ’) 


and h(/,,4,,..,@,) is the lower confidence limit for the reliability function 





h(A,,. 4, -.- a) where d is the upper confidence limit for 4,,. In his thesis, Lee [Ref. 
10: p.7] has demonstrated that if the same life data used to obtain A is used to construct 


estimates 4, for the a, in equation (1.4), then the approximate lower confidence limit 


HCA ys Bys Age oes » Ay) ~ (1.5) 


iS an accurate approximate interval estimation procedure. 

















A primary reason for investigating the accuracy of this failure rate ratio method for 
cyclic systems is that it can be easily extended to systems that have mixtures of cyclic 
and continuous components. If the procedure investigated here for systems of cyclic 

F components is found to be accurate, then it can readily be combined with the procedures 
established by Lee [Ref. 10: pp.3-23] to obtain an accurate interval estimation procedure 
for the reliability of a coherent system with a mixture of cyclic and continuous compo- 
nents. 

An additional primary purpose of this thesis is to develop an extensive computer 
simulation program that will provide a means for evaluating the accuracy of proposed 
interval estimation method. for the reliability of coherent systems with cyclic compo- 
nents. 


vag Call eid one 











Hl. INTERVAL ESTIMATION PROCEDURES FOR SYSTEM 
RELIABILITY 


A. GENERAL CONCEPTS 

Consider a coherent system with k cyclic components. Mission tests for component 
‘ consist of independent Bernoulli trials administered under a mission operating envi- 
ronment with probability g, of failure on each test. Then F,, the number of failures in 7, 
independent mission tests, has a binomial distribution denoted by BIN(n,q). We shall 
refer to gq, as the unreliability of component 7. 

Utilizing the concept first suggested by Lomnicki [Ref. 8: p. 109] and expanded by 
Myhre and others [Ref. 9: p.213], the unreliability of each component can be written as 
a fraction of the unreliability of the least reliable component; that is 


91 = 49m (2.1) 


where 


— Mm=max{g}, f=1,2,..0k and “O0<a,<1. Wé"shall assume that the q, are suf- 
ficiently small so that the probability distribution of F, can be approximated by the 
Poisson distribution with mean 7, q,. We denote this by the expression 


F,~ P(m, 9). (2.2) 


Therefore, the distribution of the total number of failures, F= oF, can be approximated 


by the Poisson distribution with mean =n, q;; that is, 
i=0 


k 
Pe (dn a) (2.3) 


i=} 
Applying equation (2.1) to equation (2.3) we obtain 
k 
jae H( ay «) (2.4) 
. d=] 


--Let h(4, %,--,4) denote the equation for system reliability, R;. Then from 


equation (2.1), 





aya 














Rg = A(Giny Qj, Ay, «++ y Ay) (2.5) 


Since the system is coherent, h is non-increasing in g,- Therefore, a lower 
100(1-a)% confidence limit, Rs 7 for Rs is 


Rs, 142 = (Gin, pays 29 «> + Me) (2.6) 


where Gnu) is any 100(1-«)% upper confidence limit for g,. To obtain Gnas WE use the 
following well known result for the upper confidence limit of the mean of a Poisson 
distribution, see [Ref. 11: p.218]. If X ~ P(Z) then the upper 100(1-«)% confidence limit 
Awa) for A using one observation on_X is 


2 
; Xa, 214%) 
ducay = — ; (2.7) 


Applying this result to equation (2.4) we obtain the upper 100(1-x)% confidence limit 
A 
for q,>.", a, Which is given by 
i=t 


k 2 
Xn, 21+F) 
@x +a Se Soa (2.8) 
i=] 
The corresponding upper confidence limit for q,, is 
7 
a, 21+F) 
Ym u(x) = k : (2.9) 
2 Ny G; 


i=] 


Equation (2.9) is the expression used by Myhre and others [Ref. 9: p.214 J. ; 

The components of the vector g = (a,, @,...,@,) aS defined by Mvhre and others 
[Ref. 9: p.214] are computed in relation to the largest unreliability g,. However any 
component could be used as the base component to form a vector of ratios. Choosing 
the one with largest unreliability is a convenient way to select the index m. In this thesis 
we will assume the gq, are unknown. We will estimate them from the data. We also use 
the data to determine the base component, which we label component m, that will be 
used to form the ratio estimates a@,. Specifically if N,, F, are the number of trials and 
resultant failures for component i, then g, = + is the estimate of unreliability ¢, Here- 
after in this thesis, we define the index m bv the equation 





phate 


iat 











am = max(q), 9a, song 9,): (2.10) 
We shall define a, by 
fee Fh ck (2.11) 
Im 


Note that 4, is not intended to be an estimate of max(q,, ..., q,). It is nothing more than 
the largest of the observed estimates of the component unreliabilities and provides us 
With the basecomponent which we denote by m. Hereafter, in this thesis, g, will denote. 


the unreliability of the component m which has been determined by the definition of g,._ 


With the index m determined and the quantities 4,, @,..., 4, computed, we can com- 
pute an approximate upper confidence limit 4,44) for g, that corresponds tO gp.) in 
equation (2.9). Specifically, 


2 
Xa, U+F) 
Gmnaute) = + : (2.12) 
A 
2 el a; 
i=} 


The corresponding approximate lower 100(1-a)%> confidence limit Rs.rey iS 
HG atayy B19 «+ » Ay): (2.13) 


We have special problems with equation (2.12) when none of the components fail. In 
this case g,=0, (i=1,2,...,4) 9, =0 and all a, are undefined. If at least one compo- 
nent: has at least one failure, equation (2.12) is well defined. Consequently, equations 


(2.12) and (2.13) must be modified or supplemented to account for cases when none or © 


perhaps only one of the components fail. In fact we may want to modify the definition 
of a, whenever @,=0; i.e. when F,=0. Anytime a,=0, the N, tests for component i 

make no contribution toward the evaluation of g,, 4) in equation (2.12). If information 
about the range of qa is known, it could be used to redefine a, so that @,n, in the 
denominator of equation (2.12) would not be zero. This is mathematically equivalent 
to adding a,n, tests to the N,, tests for component m without adding any failures to F,. 
This is the one advantage-of the ratio method and itis reflected in both equations (2.12). 
“and (2.13). It is this specific property of the ratio method that makes it particularly 
appealing when the total number of failures is small. 





rg 











The primary purpose of this thesis is to initiate the development of alternative 
methods for constructing approximate lower confidence limits Rs 1) for the reliability 
Ry of coherent systems that account for zero or nearly zero failures and to construct a 
computer program that can be used to evaluate the accuracy of these confidence limit 
procedures. It is highly desirable that these methods make strong use of the ratio 
method, when, several types.of components experiencg failure, because doing so should 
allow us to extend these confidence limit procedures to systems that have both cyclic and 
continuously operating types of components. . 

It is highly unlikely that any one confidence limit procedure will be reasonably ac- 
curate for all system configurations; i.e., for all reliability functions, h(y, ...,9,). Con- 
sequently, we have developed three similar procedures which are labeled Procedure 1, 2, 
and 3. Each of these procedures will be evaluated for accuracy when applied to two 
radically different types of system configurations; namely series system and Wheatstone 
bridge systems. Series systems have no redundant components. Wheatstone bridges are 
highly redundant systems. These accuracy evaluations should help establish some pre- 
liminary boundary constraints on the application of these procedures and provide in- 
sights for modifications that may vield more accurate procedures. 

The first evaluations will be performed for a series system of k independent compo- 


nents. In this case the system reliability R, is given by R, = rd —-4)= nu — AGn) 
where q, denotes the unreliability of component 7, and m will be the index established 
by the data as previously discussed, and a, = a . One expression that can be used to 
construct a lower 100(1-«)% confidence limit Rs 1.) for Rs usINg 9p, ya) BiVen in equation 


(2.12) is 


k ‘ 
Ria) = [ Ja — 4 Gn syay): (2.14) 
i=l 


In the following descriptions, m, and F, denote the number of tests (Bernoulli trials) 
and the number of failures, respectively, for component i,i=1,2,...,k. 


"B. PROCEDURE 1 


This procedure has three expressions for Rs. ts) Which depend on the test results 
F, F,,...,F, fno components fail, all g,=0, and all a, are undefined. This precludes 
the use of equation (2.12) to construct Rs. toy for this case. When no components fail, 
we define 











n* = min(m, my, .. Mm) (2.15) 


and interpret the data as n successful system tests. Rather than use the standard ap- 
— . es * 
proximate binomial confidence limit (a);7, we choose to define Rs. ray BY 


2 
Xx, 2 


2n* 





Rs wy =1- (2.16) 
This expression reflects the Poisson approximation to the Binomial distribution and will 
be slightly conservative when all F,=0, but the other components to Procedure | are 
suspected to be slightly optimistic. So the entire procedure may be nearly exact. If the 
entire procedure is conservative, we can always change equation (2.16). 

The second expression in Procedure 1 addresses the case when all but one compo- 
nent, say component 1], have zero failures. For a series system, this would correspond 
to at most one system failure out of n™ system tests. This would not be the casc in other 
systems if the failed component is redundant in the system. Here again we choose 
slightly conservative procedure for the lower confidence limit; namely, we define q,,.) 
by 


7 


2 
Xa, WI+F) 
A ’ 1 
4 u(a) = 2 . (2.17) 


We use equation (2.13) to define Rg ,,,, for a general coherent system. If the system is 
a series system we use equation (2.14). 
The third expression in Procedure 1 addresses the case when at least two different 


tvpes of components have failures; that is, F,#0 for at least two values of i 


9 


a 


aaa? i ae of A a 
i=1,2,...,k.. In this case we use equation (2.12) to define Ym.uey Where Pa em 
a t A A A . . . . q, 

4. => and q, = max(q, ...,9,). The lower confidence limit for the system is given by 
equation (2.13) for general coherent systems or equation (2.14) for a series system. 


C. PROCEDURE 2 

Procedure 2 differs from Procedure 1 in two respects. First, we assume that at least 
two different components of the system have at least one failure. Operationally. this 
means that Procedure 2 would only be applied to data for which at least two types of 


system components experienced one or more failures. In this procedure we always use 














ail Sal a neta 


equation (2.1) to obtain g,, w) and equation (2.13) or (2.14) to compute Rs, La) fOr gencral 
coherent systems or series systems respectively. Again Rae. but we will define g, 
differently for components with zero failures. Specifically, we apply a common scaling 
factor K to increase the sample size from n, to Kn, for all system components i which 
experienced zero failures and record 1 failure for each of those components. For ex- 
ample, if m,, m, and n, tests were performed on component types 1, 2 and 3 and no failures 
occurred on any of these tests, the data is modified to show F,= 1, fF; =1, F,=1 and 
new sample sizes Kn,, Kn, and Kn,. The scaling factor is an input parameter to the 
computer programs used to evaluate this procedure; ie., it is assumed to be known. It 
is determined by 


pe max(q), --- » Yx) _ 4max 


Min(g;,.-5+9,)  Fmin 218) 


where the q, are the component unreliabilities. Operationally, this would mean that the 
user of this procedure must determine an estimate for this ratio. Previous data is often 
available to provide this estimate. Such estimates can sometimes be constructed from 
Department of Defense (DOD) documents that provide a variety of quality data char- 


acteristics for hardware purchased in accordance with prescribed standards and specifi- 


cations; c.g., MILHDBK-217E and its referenced documents. Also the reliability values . 


for components used in reliability prediction analysis are usually derived from some of- 

ficial source Which could bé"satisfactoty for such estimates. A reliability prediction 
analvsis is usually required in any major DOD system acquisition program. Ifa Satis- 
factory estimate for K cannot be obtained, then Procedure 2 should not be used. Of 
course onc of the purposes of this thesis is to investigate the accuracy of this procedure 
when the scaling factor K is chosen “correctly”; i.c., according to equation (2.18). 
Thinking of ¢,., and g,i, aS failure rates, zt and a denote the expected number of 
tests to first failure for their respective components. Then the value K in equation (2.18) 
is the ratio of these two expected number of tests to obtain one failure for each compo- 


nent. That is, 


] 


min 


Viera 


max 


K=( 





). (2.19) 

















For each system component / that experiences zero failures, the corresponding esti- 
‘mate q, for q, is given by 





4 = Rn (2.20) 
Then 
ga (2.21) 
Ym 
and 
2 7 

A -a, 2(1+F) 
Fmula) ~~ (2.22) 

2) 14 

i=0 


where F is the number of actual failures observed before the data is adjusted. Thus the 
data is adjusted to obtain estimated values a, different from zero for those components 
that have no failures. The value for Rs, ta) IS given by equation (2.14). Note that if all 
k types of components have at least one . failure, the qe sulting value of Rs, Ly iS the same 
under Procedures 1 and 2. a 


D. PROCEDURE 3 

Procedure 3 differs from Procedure | in only one respect. We scale the sample sizes 
as we did in Procedure 2. That is, for each system component / that has no failures, we 
estimate gq, from equation (2.20). Equations for @,, Gn, and Rs Lay are given by equations 
(2.21), (2.22) and (2.13) or (2.14) respectively. 

Procedure 3 does not require that at least one failure have occurred on at least two 
different types of components as required in Procedure 2. Operationally, Procedure 3 


can be used for all sets of data including those sets where no failures occur. 


E. BRIDGE SYSTEM 
We define a Wheatstone bridge system by the reliability block diagram shown in 
Figure 1 on page 11. If p, denotes the reliability of component i, i= 1,2,..,5, and 


l—p,=4;=4@,q, Where ¢, = Max{q,, 93... , 9s}, then system reliability R, is given by 
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——©) 


Figure 1. Block Diagram of a Wheatstone Bridge Structure. 


Rs = (Gin @yy + 5 Gs) 
2 3 
= 1 = Gn(142 + 84s) — Fmn( 414395 + 220344) 
4 (2.23) 
+ Ypy( A, Q7034Q5 + A) Q,A,A5 + A) A,AgA, + A,A3A,a5 + 003;04a5) 


5 
— 24 jn, 42434445) 


(see Mvhre and others [Ref. 9: p.215])... Then the equation for the lower 100(1-a)% 
confidence limit on this system is 


Rs, L(x) = (Gin, way’ ay, sony as) (2.24) 


where q,, ,, and 4@, are defined as in Procedures 1, 2 or 3. Only Procedure 1 was evalu- 
ated for the bridge system and reported in this thesis. 








II. SIMULATION PROCEDURE 


Standard computer simulation methods were used to evaluate the accuracy of the 
three interval estimation methods discussed in section If. The input parameters needed 
to run the computer programs for each of the three interval estimation procedures are 
as follows: 

¢ k: number of components in the system 
© n=(n, 7, ...,n,): vector of sample sizes(mission tests) 
© G= (Gy 4, +5 9,): Vector of component unreliabilities 


—@ a: level of Confidence ~~ - sd A. 


The NON-IMSL uniform random number generator SRND was used to simulate 
the outcomes of mission tests on the A components. A total of >, uniform random 
numbers were generated and transformed into ones or zeroes as follows: the n, numbers 
Xiis Xiay ++» Xjq, in the ith block of uniform random numbers were transformed by the ex- 
pression 


Y,=1 if xySq 
Y,=0 if xy> 


for j=1,2,...,%, and /=1,2,...,4, where q, is the input parameter denoting the unre- 
liability of component i. Then F, = ¥{Y,, denotes the number of failures in m, independent 
Bernoulli trials. This sct of data is used to compute the values of g,, &, Ym,us) and Rs. ae 
- for each of the three interval estimation methods. Thus one realization of the lower 
bound estimate Rue for the random variable Rs. ty) 1S Obtained for each set of on, uni- 
form random numbers gencrated. Each realization described was replicated 1000 times : 
in order to generate a simulated empirical population of this random variable Rs. La) ° 
The vector of replications thus obtained was sorted from smallest to largest to obtain 
Riis Stile Riis Then the 100(1-a)% approximate lower confidence bound on the 
true system reliability R, is the 1000(1-x)th element of the sorted vector of replications; 
namely, Risse . . 

The simulated true confidence level is calculated by finding the element of the vector 


at 


of replications which is closest to R, and observing its index number j.. Then the simu- 


12 


- 




















lated “true” confidence level is Oy x 100. If there is a sequence of indices with the 
‘Same valuc closest to Rs , we chose the smallest index j. This may yield an arbitrarily 
small value for our recorded “true” confidence level. 

In addition to the external subroutines programmed by the author and to the ones 
mentioned above, the following subroutines were used: | 


¢ IMSL subroutine MDCHI was used to calculate the value of a chi-square random 
variable given the quantile and the degrees of freedom. 


° NON-IMSL subroutine SHSORT was used to sort, the arrays of system reliability 
, estimates in ascending order in order to obtain the appropriate order statistic. 


e IMSL subroutine USMNMX was used to extract the minimum and maximum 
values of the vector of mission tests for each component and of the vector of input 
component unreliabilities. 


Each set of input paramcters defines a simulation run, or a case. The following ta- 


bles specify the sets of parameters uscd for cach case that was simulated. 


Table 1. INPUT NUMBER OF MISSION TRIALS N, AND UNRELIABILITIES 
Q, FOR A SERIES SYSTEM. CASES 1-3. 


n, by Case Number . 
Input q; 


(cases 1-3) 


a a 


Component no. 


00100 
a CT 
a 





13 








Table 2, INPUT NUMBER OF MISSION TRIALS N, AND UNRELIABILITIES 
Q, FOR A SERIES SYSTEM. CASES 4-6. 






n, by Case Number 
Input q, 
_ Meese 4-6) 


Table 3. INPUT NUMBER OF MISSION TRIALS N, AND UNRELIABILITIES 
Q, FOR A SERIES SYSTEM. CASES 7-9. 







n, by Case Number 


Component no. La 
7 


710050500 
SS A Se re 









Input q; 
(cases 7-9) 









ee ed 
A COA! A NST A ATC 
i 


fee ieee ee 0.0050 | 


rea 











Table 4. INPUT NUMBER OF MISSION TRIALS N, AND UNRELIABILITIES 
Q, FOR A SERIES SYSTEM. CASES ‘10-12 gare 


n, by Case Number 
Input q, 


0.0025 





Table 5. INPUT NUMBER OF MISSION TRIALS N, AND UNRELIABILITIES 
Q, FOR A WHEATSTONE BRIDGE. CASES 13-15. 


Input g; 
(cases 13-15) 

















Table 6. INPUT NUMBER OF MISSION TRIALS N, AND UNRELIABILITIES 
Q, FOR A WHEATSTONE BRIDGE. CASES 16-18. 












n, by Case Number 


— ee 


ee aR Ea, a ee 
Ey aa a a 
pst 200.0350 
(eae Sa ES ee a (ee 
ee a a ee a el ee Se 





Input gq, 
(cases 16-18) 


The values of the scaling factor K used in Procedures 2 and 3 are as follows: 


Table 7. VALUES OF SCALING FACTOR K. 


a ee 





ae valucs are computed from the assigned input parameters for these cases; ¢. B. teas 


OS 
1 Ogg! 
4= Or from table 4 and 1.43~ 035 from table 6. 


CHA 





IV. SIMULATION RESULTS 


The results of the simulation runs for the twelve series system cases are presented in 
Tables 8 through 19. The results of the simulation runs for the six bridge system cases 
are presented in Tables 20 and 21. The accuracy of a particular interval estimation 
procedure for each case can be assessed from these tables. If the interval estimation 
procedure is exact, the two values R, and Rikeues would be equal and the True Con- 
fidence Level value would be 100(1-x). For example, consider Table 8, case 2 and 
a =.20. True system reliability, R, is 0.95572 and Rovecteng = (0.96398. Thus the sim- 
ulations indicate that the 80th percentile point of the distribution of Rays for Procedure 
2 for case 2 is 0.964 instead of 0.956. This error reflects the inaccuracy of the approxi- 
mate 80% lower confidence limit Rs. tem as defined under Procedure 2. The table also 
indicates that R, = 0.95572 is approximately the 77.4 percentile point of the distribution 
of Rs, ta instead of the 80th percentile point. 

The number of components in the system that had zero failures in each of the 1000 
replications was recorded. The average of the 1000 values is displayed in the last column 
of the table. 


wer ett Daal ad 


. A. SERIES SYSTEM 


Table 8. SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
2, CASES 1-3 





_ | Avg. no. 
No. of pie Son of comp. 


rane idence A 
Compnts - d L(a), 10001~«) phe with no 
failure 


0.95572 





0 
0.98572 100.00 














increase. 


Table 9. SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
2, CASES 4-6 


























No. of F True Con- en 
Compnts Rye), 1001-2) fidence with no 

P Level 
failure 





| a=20 | 0.97768 | 0.96737 _| 100.00 | 4.15 __| 
4.15 
FR arr cc ee 


Table 10,, SUMMARY OF RESULTS FOR A SERIES SYSTEM- PROCEDURE 
2, CASES 7-9 rae 


wae Avg. no. 

No. of : A of comp. 
Compnts Ry, 1008-0 : with no 
failure 


as ae one ae 


Le cee ome inne mee one 
Tae .20 | 093206 | 094925 | 70.00 | 9.01 
0.93206 [0.390% 
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increase. 


Table 11. SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
2, CASES 10-12 




















Avg. no. 
of comp. 
with no 











No. of 
Compnts: . 


Case x 
Rue, 100X(1—e) 








failure 


no 
so jZ=20 { osess_{ oseres | 7530] 9.03 
| a=.0s | 0.96552 | 0.94091 [100.00 | 9.03 | 
| x= .20 | 0.96552 | 0.94046 [100.00 [9.36 
12 | %=20_| 0.96552 | 0.94925 | 100.00 [9.48 | 
| a=.05 | 0.96552 | 0.90729 | 100.00 | 9.48 | 


“mvt ete tee 





Table 12. SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
































No. of . True Con- sr Ree 
Compnts Rijs, 1001-2) oe with no 











failure 

, a= 20 397 
° az=.0 
a 











0.95572 














increase. 


Table 13... SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
1, CASES 4-6 



















Avg. no. 
of comp. 
with no 





No. of Case 


Compnts Rye), 0001-«) 





failure 


no 
a pee [esr Temes [ioras aa 
a= .0s | 097768 [095954 | 10000 | ala 
Ps Peet eset assour {roo as 





=.05 
















Table 14. SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
1, CASES 7-9 














Avg. no. 
of comp. 







True Con- 
i e : 
fidenc with no 


Level failure 


eo | oos06 | oss | sa | 3as 
0.95254 
0.93206 | 0.88021 
084181 







No. of 


A 
Compnts Rie, 1000-0) 





“1 
2 





increase,” 


Table 18. SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
1, CASES 10-12 


































No. of Signif- ; True Con- nee HO 

Coin icance R fidence se oa 
pnts Level L(a), 1001-2) Vecel with no 
failure 


0.96552 | 0.84181 





Table 16. SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
3, CASES 1-3 


No. of ‘as Signif- : True Con- | AY8: 20. 


: of comp. 
: c : 
Compnts . oe Rea), 1000-2) yeh with no 


failure 


0.95572_| 0.93563 








eat 
ae 





increase, 


Table 17, SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
3, CASES 4-6 

































_ | Avg. no. 
No. of A ee ses of comp. 
Compnts R iexricoos—e) oe with no 


Level -: 





failure 


ea20 | oss | osen7 | 10000 [414 
a= .05 | 0.97768 | 0.93981 | 100.00 | 414 
100.00 | 453 
453 


Table 18. SUMMARY OF RESULTS FOR A SERIES SYSTEM-PROCEDURE 
3, CASES 7-9 



































Avg. no. 






































Signif- True Con- 
cack lance Runnin | fidence | oF con 

Level Level failure 
a= .20 
2=.20 | 0.93206 | 0.88770 
a=.05 | 0.93206 
a = .20 0.94925 
z= 08 [0.93206 





to 
tu 





Cid 











increase. 


Table 19. SUMMARY OF RESUETS FOR A SERIES SYSTEM-PROCEDURE 
. 3, CASES 10-12 


tail 
































_ | Avg. no. 
No. of A ey Con of comp. 
Compnts Rie), so001-«) es with no 


Level failure 


————————e—e 
| a=.05 | 0.96552 | 0.94091 | 100.00 | 9.02 
SS Eee 

| «=.05 | 0.96552_| 0.89162 | 100.00 | 9.35__| 


B. BRIDGE SYSTEM 


Table 20. SUMMARY OF RESULTS FOR A BRIDGE SYSTEM-PROCEDURE 
1, CASES 13-15 ae 



















































oer Avg. no. 

: Signif- True Con- 
ozo! icance R fidence | of comp. 
Compnts Level Lia), 10001~2) Level with no 










failure 


za .0S | 0.95782 
22.05 | 0987 

ee 1.37 
3.00000 
































increase. i 


Table 21. SUMMARY OF RESULTS FOR A BRIDGE SYSTEM- PROCEDURE 
1, CASES 16-18 








































nt Avg. no. 
Signif- True Con- 

No. of icance R fidence of comp. 

Compnts Level L(a), 1001-2) Level with no 















failure 


aa.20 | onus | 096122 | sao | 09a 
a= .05 [04816 | 0.56740 | 82.80 | 0.94 | 










0.98816 | 0.97125 





“viet el ad aptle 





ps2 | osasie | ossess | 79 | on : 











V. SUMMARY AND CONCLUSIONS 


The results in Tables 8 through 19 in section IV are quite mixed with respect to the 
accuracy across all twelve cases. However, these results definitely indicate that any one 
of the three interval estimation procedures offers excellent potential for a reasonably 
accurate lower confidence interval estimation method for system reliability. The accu- 
racy of approximate discrete interval estimation methods for system reliability in general 
depend | on. the amount of testing relative to the expected number of failures, En VW 
Woods ‘and Borsting [Ref. 7] for example, show thcir lower confidence limit procedure 
to be quite accurate for cases for which the expected number of failures is greater than 


5. The expected number of failures for some of the cases simulated in this thesis are as 
follows: 


Table 22. EXPECTED NUMBER OF FAILURES 


Consequently, any cases other than those numbered 13, 14, and 16 should not be ex- 
pected to be highly accurate if these three procedures have features similar to those of 
the Woods-Borsting procedure (Ref. 7: p.1]. Actually, in light of the small values of 













cae No.: 


: aM, q, for cases 1, 4, 7, and 10, it is somewhat surprising that the procedures are as ac- 
curate as they are at the 80% level of confidence. The results for cases 13, 14, and 16 
are quite good at the 80% level of confidence for the bridge system. Some of the pa- 
rameter sets for the cases run were designed to detect instances where the methods might 
not be accurate. 

The three procedures analyzed in this thesis should be modified and combined to 
establish one general procedure. The corresponding range of parameter scts for which 
the procedure js accurate should be established. This range will likely be a function of 


a 
Sng, This kind of a procedure could be useful in many practical settings. Morcover, 
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it could be combined with the method developed in the thesis by Lee (Ref. 10: pp.3-23]. 
Additional analvsis, and simulation need to be performed to establish such a procedure. 











APPENDIX A. INTERVAL ESTIMATION PROCEDURE 1 


PROGRAM ZFYSCN 
Ihididiininiiecicidncicccciccicoiccicccciciccceccccicicicccicicicccdekticicscicicteiciciciicrercicctctetescn 


* = * "TITLE: BINOMIAL INTERVAL ESTIMATION PROCEDURE 
is ZERO FAILURES ALLOWED; NO SCALING 

¥ AUTHOR: E. F. BELLINI, LT, USN 

v DATE: NOV 89 


* THIS PROGRAM COMPUTES THE TRUE CONFIDENCE LEVEL FOR THE ESTIMATE 
a RELIABiLITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY 


: OF THEIR COMPONENTS 

ve IN ITS PRESENT CONFIGURATION THIS PROGRAM IS SET UP TO RUN 12 

¥ TIMES EACH TIME PRODUCING 1000 REPLICATIONS USING A DIFFERENT 

¥ SET OF INPUT DATA. RUN THE PROGRAM FROM CMS BY TYPING 'B1 EXEC’ 

* THE REXX EXEC PROGRAM 

¥e 'B1' CALLS THE INPUT FILES“TO BE READ AND NAMES THE. 12 OUTPUT ove 


ve FILES RESULTING FROM THE 12 CONSECUTIVE RUNS. BY EDITING THE 
INDEX COUNTERS I, J, K OF THE ‘Bi’ ©XEC ONE CAN RUN ANY USER- 
SPECIFIC RUN FROM JUST ONE RUN TO ALL 12. 


LESSER SERRE EE SEES 


VARIABLES USED 


AHATI : WEIGHT ESTIMATES FOR EACH COMPONENT * 
. AT : INPUT WEIGHTS FOR EACH COMPONENT Big! 
ALFA  : LEVELS OF SIGNIFICANCE ¥ 
BIGF : TOTAL NO. OF FAILURES FOR EACH REPLICATION ¥ 
CHISQ : CHI-SQUARE RANDOM VARIABLE VALUE * 
C1C15 : FORMAT LABEL * 
DEGFR : DEGREES OF FREEDOM * 
DELBRG : DIFFERENCE FOR BRIDGE SYSTEM ¥e 
DELSTR : DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM ¥ 
DELTAR : DIFFERENCE FOR SERIES SYSTEM ¥ 
DIFF : DIFFERENCE (TRUE REL. - ESTIMATED REL. ) ¥ 
EPS : SMALL QUANTITY( CONSTANT) ¥ 
ERROR : PARAMETER FOR IMSL ROUTINE ¥ 
FAILS : COUNTS NO. OF REPLICATIONS WITH AT LST. 1 FAILURE * 
PI : NO. OF FAILURES FOR EACH COMPONENT(ALL MISSION TST)* 
FLAG : 1 IF ALL COMP. HAVE SAME NO. OF MISSION TESTS * 
INC : INCREMENT STEP SIZE FOR ROUTINE USMNMX * 
_ KEY1 =: ARRAY OF INDECES FOR ROUTINE SHSORT _ * 
‘KEY2  : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KEY3 : ARRAY OF INDECES FOR ROUTINE SHSORT 
KEYS : ARRAY OF INDECES FOR ROUTINE SHSORT ¥ 
KK : ARRAY SIZING PARAMETER FOR THE MAX NO OF Ce ENTS 
LOOP : COUNTS NO. OF REPLICATION PERFORMED 

MAXALF : MAX NO. OF SIGNIFICANT LEVELS DESIRED( ARRAY SIZING) 
MAXREP : MAX NO. OF REPLICATIONS 

MAXRUN : MAX NO. OF PROGRAM ITERATIONS ALLOWED 

MSTRQ : MASTER UNRELIABILITY(USED WITH AI'S TO CALC. QI'S) * 
MULT : MULTIPLIER FOR RANDOM NO. GENERAPOR SRN: 














* N : NO. OF MISSION TEST FOR EACH COMPONENT * 


* = sNIMAX =: MAX NO. OF MISSION TESTS * 
a ey ~“NIMIM : MIN NO. OF MISSION TESTS Senn 
* NINDX : INDEX NO. OF MAX NO. OF MISSION TESTS ~ 
* NIREAL : NO. OF MISSION TESTS TRANSFORMED TO REAL * 
* NMAX : MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL * 
. NPRNT : FLAG FOR DETAILED REPORT OUTPUT * 
* PRNT : SAME AS ABOVE( PARAMETER) * 
= QHATI : UNRELIABILITY ESTIMATES FOR EACH COMPONENT ¥ 
* QHTMAX : LARGEST QHATI * 
as QHTUPR : UPPER LIMIT ON SYSTEM UNRELIABILITY * 
bid QI : INPUT UNRELIABILIY FOR EACH COMPONENT * 
x QINDX : INDEX * 
* QUANTL : QUANTILE * 
* REPSHD : REPLICATIONS HEADING FORMAT NUMBER * 
* “  ““RHTSTR : SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) Fen 
bad RS =—s::~ TRUE SERIES SYSTEM RELIABILITY * 
% RSBRDG : TRUE BRIDGE SYSTEM RELIABILITY * 
* RSHAT : SERIES SYSTEM RELIABILITY ESTIMATE * 
* RSHTBR : BRIDGE SYSTEM RELIABILITY ESTIMATE ¥ 
* SEED : PARAMETER * 
* SELCTA : SIGNIFICANCE LEVEL SELECTION * 
| ve SELCTB : SIGNIFICANCE LEVEL SELECTION ¥ 
| * SORT : PARAMETER FOR ROUTINE SRND * 
* SUMNAI : SUM OF THE PRODUCT OF NI'S AND AI'S iad 
* TEMP : TEMPORARY ARRAY * 
am TOTREP : TOTAJ. NUMBER OF PROGRAM ITERATIONS * 
* TRANBR : TEMPORARY ARRAY * 
* TRANSQ : TEMPORARY ARRAY * 
¥ TRANSR : TEMPORARY ARRAY * 
al TRIALS : BERNOULLI TRIALS ARRAY (2-DIM) “ 
bad TRNSTR : TEMPORARY ARRAY * 
ba) TRUQNT : TRUE QUANTILE - 
a UNIRV : UNIFORM RANDOM DEVIATES (2-DIM) * 
* ZFAILS :- TOTAL NUMBER OF REPLICATIONS WITH ZERO FAILURES * 
¥ ZFPREP : NO. OF COMPNTS. WITH ZERO FAILURES PER REPLICATION * 
ve * 


Letedededtedetetededestede dete dedecede de sete levee ve Vee TCR TVR IRV ENEDE RE LCL NE VEDA DACA ANAS ATR AD ND LD SNA IS SNARES 


_ PARAMETER (KK=10 ,MAXALF=2 ,NPRNT=0) 

“PARAMETER (MAXREP=1000, MAXRUN=2000, EPS=. 000001) 
REAL*4 UNIRV(15,1000) , TEMP( 1000) ,QI(KK) , AICKK) , AHATI(KK) 
REAL*4 QHATI(KK), NMAX, NNMAX, QHTMAX, CHISQR(5,5), ALFA(MAXALF) 
REAL*4 DF(5),AALFA(5) , SUMNAI , RSHAT(MAXALF , MAXREP) , RS 
REAL*4 KEY1(MAXREP) ,KEY2(MAXREP) , KEY3(MAXREP) , TRNSTR(MAXREP ) 
REAL*4 DEGFR(MAXREP), QHTUPR(MAXALF ,MAXREP) , CHISQ(MAXALF , MAXREP) 
REAL*4 QUPA1(MAXREP), QUPA2(MAXREP) , RETSTR(MAXALF , MAXREP) 
REAL*4 DELTAR(MAXALF), TRANSQ(MAXREP) , TRANSR(MAXREP) , DIFF(MAXREP) 
REAL*4 DELSTR(MAXALF) ,NIMIN,NIMAX , NIREAL(KK) 
REAL*4 RSHTBR(MAXALF ,MAXREP) , DELBRG(MAXALF) , XEY4(MAXREP) 
REAL*4 TRANBR(MAXREP), RSBRDG ,MSTRQ 
REAL*4 ZFPREP 


“INTEGER SEED, MULT, SORT, TRIALS(15,1000), BIGF, FICKK), NCKK) .°™ 
INTEGER NINDX, QINDX, ERROR, REPS, SELCTA, SELCTK, TOTREP 
INTEGER C1C15, REPSHD, SELCTB, ALF, FLAG, LOOP,PRNT 

















INTEGER QUANTL(MAXALF), TRUQNT(MAXALF) ,ZFAILS ,FAILS , INC 


DATA SEED/123457/, MULT/1/, INC/1/ 

DATA AALFA/.01,.05,.9,.95,.99/, DF/1,5,10,30,40/ 
DATA ALFA/. 20,.050/ 

DATA SORT/0/ 


ASSIGN 8 TO C1C15 
ASSIGN 9 TO REPSHD 


* CALL COMPRS | a as — &. we 
PRNT = NPRNT 


DO 12 I=1,KK 
AI(I) = 9999. 
NCI) = 99999999 
12 CONTINUE 


READ(03,*)K,MSTRQ 


DO 11 I=1,K 
READ(03,*) AICI) ,NCI) 
11 CONTINUE 


IF(K.NE.5) THEN 
WRITE(1,'(''WARNING: BRIDGE STRUCTURE '', 
+''ONLY USES THE FIRST 5 COMPONENTS'')') 


ELSE 

END IF 
wk / / INITIALIZE THE QHTUPR ARRAY OF UNRELIABILITY REPLICATIONS, = //*** 
* RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS ¥ 
ve AND RHTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN * 
weewe/ / ALL THE COMPONENT MISSION TESTS ARE EQUAL IN NUMBER [ [Re 


DO 172 ALF=1,MAXALF 
DO 173 REPS=1,MAXREP 

QHTUPR(ALF ,REPS) = 0 
RSHAT(ALF REPS) = 0. 
RHTSTR(ALF ,REPS) = 
RSHTBR(ALF REPS) = 

173 CONTINUE 

172 CONTINUE 


0 
0 


week / 1 SET FLAG TO 1 IF ALL COMPONENTS HAVE SAME NO. OF MISSION TESTS *%%¥ 


FLAG=1.. 
DO 50 I=1,K -1 
IF((N(I) - NC(I+1)).NE. 0) THEN 
FLAG=0 


eos uid a 


ELSE 
END IF 
50 CONTINUE 
PRINT *, 'FLAG IS:', FLAG : 
wick 1 / MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION) / /** 











“10 


weer | | 


eae ] / 


95 


eter | | 


eae | / 


Yerek ff 


115 


125 
120 


aesere ] / 


135 
130 


ZFPREP = 0. 


ZFAILS = 0 

FAILS = 0 

TOTREP = 0 im am 

LOOP = ‘ as 


IF( LOOP. LT. MAXREP) THEN 
LOOP = LOOP + 1 
IF(TOTREP. LT. MAXRUN) THEN 
TOTREP = TOTREP + 1 
SELCTA = 1 
SELCTB = 2 
FILL ARRAY KEY(REPS) WITH INTEGERS 1 TO K TO BE USED AS OUTPUT 


OF THE SUBROUTINE SHSORT 
DO 95 REPS=1, MAXREP 
KEY1(REPS) = REPS 
KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEY4(REPS) = REPS 
CONTINUE 


CALCULATE NMAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE 
CALCULATE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NO. //*** 
CALL IMAX(N,K,NMAX,NINDX) 


CALCULATE THE QI'S FROM THE GIVEN MASTER Q AND THE AI'S 


“po 115 I=1, K 


QI(1) = MSTRQ * AI(1I) 
CONTINUE 


DO 120 I=1,15 
DO 125 J=1,500 
UNIRV(I,J) = 999. 
TRIALS(I,J) = 99999 
CONTINUE 
CONTINUE 


DRAW UNIFORM (0,1) RV'S AND CONVERT TO BERNOULLI TRIALS //*** 
‘: nee wwe tel ny 
DO 130 I=1, K 
CALL SRND(SEED, TEMP, N(I), MULT, SORT) 
DO 135 J=1, NCI) 

UNIRV(1I,J) = TEMP(J) 

IF (UNIRV(I,J).LE. 1 - QI(I)) THEN 

TRIALS(I,J) = 0 


TRIALS(I,J) = 1 
END IF 


CONTINUE 
CONTINUE 
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taker ff 


152. 


Were / | 


200 


205 


141 


146 





CALCULATE THE NO: OF FAILURES FOR EACH™COMPONENT //** 


DO 150 I=1, K 
FI(I) = 
CONTINUE 


CALCULATE THE F SUB I'S AND THE GRAND TOTAL NO. OF FAILURES 
DO 155 I=1, K 
DO 160 J=1, N(I) 
FI(I) = FI(I) + TRIALS(I,J) 
CONTINUE 
IF(FI(1I).EQ.0) THEN 
ZFPREP = ZFPREP + 1 
ELSE 
END IF 
CALCULATE THE QHAT SUB I'S: F SUB I'S DIVIDED BY N SUB I'S 
QHATI(I) = REAL(FI(I)) / N(I) 
BIGF = BIGF + FI(I) 
CONTINUE 


CASE WHERE NO COMPONENTS HAVE ANY FAILURES f [eee 


IF(BIGF. EQ. 0) THEN 
_ FAILS = Z2FAILS + 1 
“DO 200 I=1, K 
NIREAL(I) = REAL(N(1)) 
CONTINUE 
CALL USMNMX(NIREAL,K, INC ,NIMIN,NIMAX) 
DO 205 ALF=1, MAXALF 
CALL MDCHI(1 - ALFA(ALF),2. ,CHISQCALF,LOOP) , ERROR) 
RSHAT(ALF ,LOOP)= 1 - (CHISQCALF,LOOP) / REAL(2 %* NIMIN)) 
IF(FLAG. EQ. 1) THEN 
RHTSTR(ALF , LOOP)=1-(CHISQ(ALF LOOP) / REAL(2 * N(1))) 
ELSE 
END IF 
CONTINUE 
IF(PRNT. EQ. 1) THEN 
-  “WRITE(1,0007) ae = 
WRITE(1,C1C15) 
WRITE(1,3334) QI 
WRITE(1,0001) 
WRITE(1,C1C15) 
DO 141 J=1,NMAX 
WRITE(1,1111) (UNIRV(I,J), I=1, K) 
CONTINUE 
WRITE(1,0002) 
WRITE(1,C1C15) 
DO 146 J=1,NMAX 
WRITE(1,2222) (TRIALS(I,J), I=1, K) 
CONTINUE 
WRITE(1,0003) 
WRITE(1,C1C15) 
WRITE(1,3333) FI 
WRITE(1,0005) 
WRITE(1,C1C15) 


3] 











Wek / / 


week | f 


140 


145 


Heiek | / 


165 





WRITE(1,3335) N 
WRITE(1,0004) 
WRITE(1,C1C15) 
WRITE(1,3334) GHATI | 

WRITE(1,'(/'!THE MAXINUM Q HAT SUB I IS:'', T40, F8.5)') QHTMAX 

WRITE(1,'(/''THE MAXI Q HAT SUB I IS ELMNT'NO.:'' T40,15)*) QINDX 


EE Ae aes GRAND TOTAL NO. OF FAILURES IS:'',T40, I5)') BIGF 


ENDIF 
DEGFR(LOOP) = 2. 
GO TO 10 
ELSE 
FAILS = FAILS + 1 
END IF 


FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES 
CALL RMAX(QHATI, K, QHTMAX, QINDX) 


PRINT THE RESULT OF THE MAIN OPERATING ELEMENTS OF THE PROGRAM 


IF(PRNT. EQ. 1) THEN 
WRITE(1,0007) 
_. WRITE(1,C1C15) 
~ WRITE(1,3334) QI 
WRITE(1,0001) 
WRITE(1,C1C15) 
DO 140 J=1,NMAX 
WRITE(1,1111) CUNIRV(I,J), I=1, K) 
CONTINUE 
WRITE(1,0002) 
WRITE(1,C1C15) 
DO 145 J=1,NMAX z : 
WRITE(1,2222) (TRIALS(I,J), I=1, K) 
CONTINUE 
WRITE(1,0003) 
WRITE(1,C1C15) 
“WRITE(1,3333) FI “™ rene 
WRITE(1,0005) 
WRITE(1,C1C15) 
WRITE(1,3335) N 
WRITE{1,0004) 
WRITE(1,C1C15) 
WRITE(1,3334) QHATI 
WRITE(1,'(/''THE MAXIMUM Q HAT SUB I IS:'', T40, F8.5)') QHTMAX 
WRITE(1,'(/''THE MAXI Q HAT SUB I IS ELMNT NO.: ' »T40,15)") QINDX 
WRITE(1,'(/''THE GRAND TOTAL NO. OF FAILURES IS:'' ,T40, 15)') BIGF 
ELSE 
ENDIF 


CALCULATE THE AHAT SUB I'S (WEIGHT ESTIMATES) 


SUMNAI = 0. 
DO 165 I=1, K 
AHATICI) = QHATICZ) / QHUTMAX 
SUMNAI = SUMNAI + NCI) * AHATICI) 
CONTINUE 
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IF(PRNT. EQ. 1) THEN 
WRITE(1,0006) 


WRITE(1,C1C15) ee 


WRITE(1,3334) AHATI 
ELSE 
END IF 


week // CALCULATE 1 REPLICATION OF UPPR ALFA C.L. ON SYSTEM RELIABILITY 
DEGFR(LOOP) = 2 * (1 + BIGF) 


DO 170 ALF=1, MAXALF 
CALL MDCHI(1 - ALFACALF) ,DEGFR(LOOP) ,CHISQ(ALF,LOOP), ERROR) 
QHTUPR(ALF,LOOP) = CHISQ(ALF,LOOP) / (2 * SUMNAT) 
IF(FLAG. EQ. 1) THEN 
es RHTSTR(ALF ,LOOP) = 1 -(CHISQ(ALF,LOOP) / REAL(2*N(1))) 

os sesEND IF 
eg (ALF, LOOP), ALFA(ALF) 
wie// CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR COMPNTS. IN SERIES 


CALL RHTSRS(QHTUPR(ALF,LOOP), AHATI,K, RSHAT(ALF, LOOP) ) 
¥ +T40,F8.5)') RSHAT(ALF , LOOP) 


wee // CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE *** 


CALL RHTBRG(QHTUPR( ALF , LOOP) , AHATI ,K, RSHTBR( ALF , LOOP) ) 
170 CONTINUE 


wee / / THIS ELSE AND ENDIF ARE FOR THE TEST AGAINST MAXRUN ittiekdeinkieink 
ELSE ~ med 
WRITE(1, van "1 /* PROGRAM EXCEEDED THE MAX NO. OF RUNS'', 
+'' ALLOWED OF: '',I6)') TOTREP 
GOTO 9999 


END IF 
GOTO 10 
END IF 
wink // SORT THE ARRAYS OF SYSTEM UNRELIABILITIES(1 FOR EACH CONF. LEVEL) 


DO 700 ALF=1, MAXALF 
DO 800 REPS=1, MAXREP 


TRANSQ(REPS) = QHTUPR(ALF,REPS) 

TRANSR(REPS) = RSHAT(ALF REPS) 

TRNSTR(REPS) = RHTSTR(ALF,REPS) 

TRANBR(REPS) = RSHTBRCALF,REPS) 
800 CONTINUE 


CALL SHSORT(TRANSQ,KEY1,MAXREP) 

CALL SHSORT(TRANSR, KEY2 ,MAXREP) 

CALL SHSORT(TRNSTR ,KEY3 ,MAXREP) 

CALL SHSORT( TRANBR,KEY4 ,MAXREP) 

DO 900 REPS=1, MAXREP 
QHTUPR(ALF,REPS) = TRANSQ(REPS) 
RSHAT( ALF ,REPS) = TRANSR(REPS) 
RHISTR(ALF,REPS) = TRNSTR(REPS) 
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RSHTBR(ALF ,REPS) = TRANBR(REPS) 
900 CONTINUE 
700 CONTINUE 


wiee/ 7? PRINT OUTPUT REPORT HEADINGS *itietekinkieiiiric 


WRITE( 1, 6666) 
WRITE(1,6667) MAXREP 
WRITE(1,6668) K 
WRITE( 1, 6669) 
. IF(K.EQ.5) THEN 
WRITE(1, 6699) 
. < “ELSE 
END IF 
WRITE(1,6670) MSTRQ 
- WRITE(1,6671) 
WRITE(1,01015) 
WRITE(1,3334) AI 
WRITE(1, 0007) 
WRITE(1,€1015) 
WRITE(1,3334) QI 
WRITE(1,0005) 
WRITE(1,C1C15) 
WRITE(1,3335) N 
WRITE(1,6674) 
WRITE(2,6666) st newt 
WRITE(2,6667) MAXREP 
WRITE(2,6668) K 
WRITE(2,6669) 
IF(K. EQ. 5) THEN 
WRITE( 1, 6699) 
ELSE 
END IF 
WRITE(2,6670) MSTRQ 
WRITE(2,6671) 
WRITE(2,C1C015) 
WRiTE(2,3334) AI 
WRITE(2,0007) 
WRITE(2,C1C15) 
WRITE(2,3334) QI 
WRITE(2,0005) 
WRITE(2,€1C15) 
WRITE(2,3335) N 
WRITE(2,6674) 


WRITE(2,'(''SORTED RSHAT 1 IS:'',/10(F8.5))') 
+(RSHAT(1,REPS), REPS=1, MAXREP) 
WRITE(2,'(''SORTED RSHAT 2 IS:'',/10(F8.5))') 
+(RSHAT(2,REPS), REPS=1, MAXREP) 
IF(FLAG. EQ. 1) THEN 
- oss. WRITE(2,'C''SORTED RHTSTR 1 1s: '',/10(F8.5))') 
+(RHTSTR(1,REPS), REPS=1, MAXREP) 
WRITE(2,'(''SORTED RHTSTR 2 IS:'',/10(F8.5))') 


+(RHTSTR(2,REPS), REPS=1, MAXREP) 
ELSE 








END IF 
IF(K.EQ. 5) THEN 
WRITE(2,'¢''SORTED RSHTBR 1 IS:'',/10(F8.5))') 
+(RSHTBR(1, REPS) , REPS=1, MAXREP) 
WRITE(2,'(''SORTED RSHTBR 2 Is:'',/10(F8.5))') 
+(RSHTBR(2,REPS), REPS=1, MAXREP) 
ELSE 
END IF 


wiet// COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL. FNCTN. (SERIES SYSTEM) 
wiek// AND FOR THE 5-COMPONENT BRIDGE STRUCTURE 
CALL RSRS(QI, K,RS) 
WRITE(1,'('' '',////''THE TRUE SERIES SYSTEM '', 
+' "RELIABILITY VALUE IS: '' ,T51,F8.5)') RS 
CALL RBRIDG(QI,K,RSBRDG) 
IF(K. EQ. 5) THEN 


WRITE(1,'('' '',////''THE TRUE BRIDGE STRUCTURE '', 
+' RELIABILITY VALUE IS:'',151,F8.5)') RSBRDG 
ELSE 
END IF 
WRITE(1,6675) 
wee // COMPUTE THE DIFFERENCE 'DELTAR' BTWN-”RS AND RSHAT OF THE THEO . 








“tele® 77 RETICAL QUANTILE GIVEN BY ALFA(MUST USE SORTED RSHAT ARRAY) 


IF(FLAG. EQ.1) THEN 
WRITE(1,5755) 


ELSE 
END IF ; 
S DO 450 ALF=1, MAXALF 
QUANTL(ALF) = MAXREP * (1 - ALFA(ALF)) 
DELTAR(ALF) = RS - RSHAT(ALF,QUANTL( ALF) ) 
DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTL( ALF) ) 


IF(FLAG. EQ. 1) THEN 
DELSTR(ALF) = RS - RHTSTR(ALF,QUANTL(ALF)) 
WRITE(1,5555) MAXREP, ALFA(ALF), REAL( QUANTL(ALF)) 
WRITE(1,5656) RHTSTR(ALF,QUANTL( ALF) ) 
WRITE(1,5657) DELSTR(ALF) 
ELSE 
END IF 
IF(K.EQ.5) THEN 
DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTL(ALF)) 
WRITE(1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF) ) 
WRITE(1,5666) RSHTBR(ALF ,QUANTL(ALF) ) 
WRITE(1,5667) DELBRG(ALF) 
ELSE 
END IF 
WRITE(1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 
WRITE(1,5556) RSHAT( ALF, QUANTL( ALF) ) 
WRITE(1,5557) DELTAR(ALF) 
450 CONTINUE 
PRINT *, 'QUANTL(1) IS:', QUANTL(1) 
PRINT *, 'QUANTL(2) IS:', QUANTL(2) 


= weve / / FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*%%* 
Lede: / / dededededeveveve RSHAT devedslelveveds / [wiek 











WRITE(1, 6676) 
DO 400 ALF=1,MAXALF 
TRUQNT(ALF) = 0 

DO 500 REPS=1, MAXREP 


_osose,, DIFF(REPS) = RS - RSHAT(ALF,REPS) 
500 CONTINUE 
DO 600 REPS=1, MAXREP 
IF(ABS(DIFF(REPS)).LE.EPS) THEN 
TRUQNT( ALF) = REPS 
WRITE(1,'('' '',/''TRUE CONFIDENCE LIMIT Is:'' 


+ F8.4)') 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
GO TO 620 


ELSEIF(DIFF(REPS). LT. 0.) THEN 
TRUQNT(ALF) = REPS 


GO TO 610 
ELSE 
owe «= CSO<EENDSOSIF OL es sis 
600 CONTINUE 
610 IF(TRUQNT( ALF). EQ. 0.) THEN 
WRITE(1,4443) ALFACALF) 
WRITE(1,'('' '',/''THE SMALLEST'', 
+ '’ DIFFERENCE BETWEEN RS AND RSHAT Is:'' ,F10.5)") DIFF( 
+ MAXREP ) 


ELSEIF( TRUQNT( ALF). EQ. 1.) THEN 
WRITE(1,4442) ALFACALF) 
WRITE(1,'C'' '',/''ALL RSHAT'', 
+ '' ARE GREATER THAN RS'')’) 
ELSEIF(ABS(DIFF( TRUQNT( ALF) )). LE. ABS(DIFF(TRUQNT(ALF) - 1))) 
- THEN 
WRITE(1,4444) ALFACALF), 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
WRITE(1,4445) RSHAT(ALF , TRUQNT( ALF) ) 
WRITE(1,4446) 
ELSE 
WRITE(1,4444) ALFACALF), 
+ ((TRUQNT(ALF)-1) / REAL(MAXREP)) * 100. 
WRITE(1,4445) RSHAT(ALF , TRUQNT( ALF) -1) 
WRITE(1,4447) 
620 END IF 
400 CONTINUE 


éék// FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*#* 
sek { / deteieteicieiek: RSHTBR (BRIDGE) wives ] [eter 


IF(K.EQ.5) THEN 
DO 401 ALF=1,MAXALF 
TRUQNT(ALF) = 0 
DO 501 REPS=1, MAXREP 
DIFF(REPS) = RSBRDG - RSHTBR(ALF,REPS) 
501 CONTINUE 
DO 601 REPS=1, MAXREP 
IF(ABS(DIFF(REPS)). LE. EPS) THEN 
: TRUQNT(ALF) = REPS 
= WRITE(i,'(+' '',/''TRUEsCONFIDENCE LIMIT Is:'', 
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a eae F8.4)') 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
GO TO 621 


ELSEIF(DIFF(REPS). LT. 0.) THEN 
TRUQNT( ALF) = REPS 
GO TO 611 
ELSE 
END IF 
601 CONTINUE 
611 IF(TRUQNT( ALF). EQ. 0.) THEN 
WRITE(1, 4443) ALFA( ALF) 
WRITE(1,'('' '',/''THE SMALLEST'', 
+ a DIFFERENCE BETWEEN RSBRDG AND RSHTBR Is:'' 
5 oad F10.5)*+) DIFF(MAXREP) Saad 
ELSEIF(TRUQNT( ALF). EQ. 1.) THEN 
WRITE(1, 4442) ALFA( ALF) 


WRITE(1,'('' '',/''ALL RSHTBR'', 
+ '' ARE GREATER THAN RSBRDG'')' } 
ELSEIF(ABS( DIFF(TRUQNT( ALF) )). LE. ABS(DIFF(TRUQNT(ALF) - 1))) 
+ THEN 
WRITE(1,4444) ALFACALF), 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 


WRITE(1,4449) RSHTBR( ALF, TRUQNT( ALF) ) 
WRITE(1,4446) 
ELSE 
WRITE(1,4444) ALFACALF), 
+ ((TRUQNT(ALF)-1) / REAL(MAXREP)) * 100. 
WRITE(1,4449) RSHTBRCALF , TRUQNT(ALF)~-1) 
WRITE(1,4447) 
621 END IF 
- 401 CONTINUE 
ELSE 
END IF 


week / / FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*#* 
aevere J / Setedevedokicre RHTSTR weitededeicdeicte | [Ri% 


~ TF(FLAG. EQ. 1) THEN 
DO 4400 ALF=1,MAXALF 
TRUQNT(ALF) = 0 
DO 5500 REPS=1, MAXREP 
DIFF(REPS) = RS - RHTSTR(ALF,REPS) 
5500 CONTINUE 
DO 6600 REPS=1, MAXREP 
IF(ABS(DIFF(REPS)). LE. EPS) THEN 
TRUQNT(ALF) = REPS 
WRITE(1,'('' '',/''TRUE CONFIDENCE LIMIT IS:'', 


+ F8.4)') 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
; GO TO 6620 
> ELSEIF(DIFF(REPS).LT.0.) THEN ; o 
TRUQNT(ALF) = REPS 
GO TO 6610 
ELSE 


END IF 














ee _ 
6600 CONTINUE Pe 
6610 IF(TRUQNT( ALF). EQ. 0.) THEN . 
WRITE(1,4443) ALFA( ALF) 
WRITE(1,'('' '',/" THE SMALLEST'', 
+ '' DIFFERENCE BETWEEN RS AND RHTSTR IS:'' 
" * F9.5)') DIFF(MAXREP) 
ELSEIF(TRUQNT( ALF). EQ. 1.) THEN 
WRITE(1,4442) ALFA(ALF) 
- owe WRITEC1y'('' 1 od" "ALL RHTSTRI', a 
+ ' ARE GREATER THAN RS'')') 
‘ ELSEIF(ABS(DIFF(TRUQNT( ALF)))- LE. ABS(DIFF(TRUQNT( ALF) -1))) 
WRITE(1,4444) ALFA(ALF), 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
WRITE(1,4448) RHTSTR(ALF , TRUQNT( ALF) ) 
WRITE(1,4446) 
ELSE 
WRITE(1,4444) ALFA(ALF), 
+ ((TRUQNT(ALF)-1) / REAL(MAXREP)> * 100. 
WRITE(1,4448) RHTSTR(ALF , TRUQNT( ALF) -1) 
WRITE(1,4447) 
| 6620 END IF 
4400 CONTINUE 
| ELSE 
END IF 


wee / / PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION *iriev 
IF( PRNT. EQ. 1) THEN 
I=1 


185 WRITE(1,REPSHD) ALFA(SELCTA), ALFA(SELCTA), 
+ALFA(SELCTB) , ALFA( SELCTB) , ALFA( SELCTA) , ALFA(SELCTA) , ALFA( SELCTB) , 
+ALFA( SELCTB) 
175 .-1FCI.GE. (MAXREP + 1)) THEN 
GOTO 180 
ELSE 
IF( (I.EQ. 71). OR. (I. EQ. 211). OR. (1. EQ. 351). OR. (1. EQ. 491). OR. 
Es (1. EQ. 631). OR, (1. EQ. 771). OR. (1. EQ. 911). OR. (1, EQ. 1051) ) THEN 
I = 1+ 70 
WRITE(1, "C04"? 9" 
GOTO 185 
ELSE 


WRITE(1,3336) I, INT(DEGFR(I)), CHISQ(1,1I), QHTUPR(1,1), 
+ CHISQ(2,1), QHTUPR(2,I) 


END IF 
“IF(I + 70. LE. MAXREP}*THEN “one ove 
WRITE(1,3337) I+70, INT(CDEGFR(I+70)) ,CHISQ(1,1+70), 
+ QHTUPR(1, I+70), CHISQ(2, I+70) ,QHTUPR(2,1+70) 
ELSE 
END IF 
T=I¢+1 
GOTO 175 
180 END IF 
ELSE 
ENDIF 














9999 WRITE(1,'(''THE TOTAL NO OF REPS WAS: '',I8)') TOTREP 
WRITE(1,'(''THE TOTAL NO OF EFFECTIVE REPS WAS: 18)’ ) LOOP 
WRITE(1,'(''THE TOTAL NO OF NO FAILURE RUNS WAS: 7 ,18)') ZFAILS 
WRITE(1, "(*' AVERAGE NO. OF COMPONENTS PER REPLICATION WITH ty 

+. "NO« FAILURES: 't4,F5.2)')}ZFPREP / MAXREP 

: WRITE(1,'' C1 THE TOTAL NO OF RUNS WITH FAILURES WAS: "',I8)') FAILS” 

0008 FORMAT ¢/ 3X,' C 1' ,5X,' C 2" 

+5X, 'C 3! 25X, C 4' 5X, Cc 5' 5X, 8 6" ,5X,'C 7' 5X, 
+'C 8' SX: "C 9! 29K, ‘Cc 10' 4X, "C 11' »4X, 
+'C 12° ,4X, C 13" 4X, "¢ 14" 4X, "C 15! ) 

0009 FORMAT( /1X, "REP No" »2X," DF 1X, "CHISQR(', F4.3,')' 2 1X, 8, 
+! QHTUPR( ' »F4. 3,! baa 1, CHISQR( re 3 ie de ee  QHTUPR( »F4.3,')' 
+2X, "REP No' s2X5 ‘DF , 1X  CHISQR( 3; oy 
+! QHTUPR(' FA, eee i ,1X, CHISQR(' fs. 3, 3K,  QHTUPR(' »F4.3,')'/) 

0001 FORMAT C7) UNIFORM RANDOM DEVIATES ine: ry 

0002 FORMAT C///' BERNOULLI TRIALS ARE: 

0003 FORMAT (///'TOTAL NO. OF FAILURES FOR EACH COMPONENT: ), 

0004 FORMAT (///'ESTIMATED UNRELIABILITY FOR EACH COMPONENT: " ) 

0005 FORMAT (///'TOTAL NUMBER OF MISSION TESTS: ') 

0006 FORMAT (///'ESTIMATED WEIGHTS FOR EACH COMPONENT: ' ) 

0007 FORMAT (///'Q I FOR EACH COMPONENT: ') 

(15F 


1111 FORMAT 8.5) 


2222 FORMAT (/1X,15(14,4X)) 

3333 FORMAT (/1X,15(14,4X)) 

3334 FORMAT (/15F8. 5) 

3335 FORMAT (/1X, Ce 4X)) 

3336 FORMAT (T3,14,T9, 13, T13,F11.5,T27,F8.5,1T39,F11.5,T53,F8. 5) 

3337 FORMAT Ce 767, 14, 173,13,T77,F1l. 5,791, F8.5, T103, F1l.5,1117,F8. 5) 

4442 FORMAT UT THE RESULTING (1 - ',F4.3,') CONFIDENCE ', 
+" LIMIT ',T50,' 00.000 ') 

4443 FORMAT Sie THE RESULTING (1 - ',F4.3,') CONFIDENCE ', 

+'LIMIT ,T50,'100. 0000" ) 
4444 FORMAT Ve THE RESULTING (1 - ',F4.3,') CONFIDENCE ', 


( 
I 
( 
I 
C 
+'LIMIT I 
4445 FORMAT ( mh THE RSHAT VALUE CLOSEST TO RS IS: ",T51,F8.5)° 
4446 FORMAT (' ',/'(FIRST NEGATIVE DIFFERENCE) ' ) 
(' ',/'CELEMENT PRECEEDING FIRST NEGATIVE DIFFERENCE )' ) 
( vol THE RHTSTR VALUE CLOSEST TO RS IS: 151, F8. 5) 
( ad THE RSHTER VALUE CLOSEST TO RSBRDG IS: »T51,F8.5) 
( pd THE ',14,'(1-',F4.3,') QUANTILE IS: ,T49, F8. 3) 
(' ',/'THE VALUE OF RSHAT FOR THAT QUANTILE IS: ,T51, F8.5) 
( ,/'THE DIFFERENCE(RS - RSHAT) IS: ,T51,F8.5) 
5656 FORMAT { st THE VALUE OF RHTSTR FOR THAT QUANTILE IS: ',T51,F8. 5) 
5666 FORMAT ( ie THE VALUE OF RSHTBR FOR THAT QUANTILE Is: ' 751, F8.5) 
5657 FORMAT (' ',/'THE DIFFERENCE(RS ~- RHTSTR) IS: ,T51,F8. 5) 
5667 FORMAT ( ie THE DIFFERENCE(RS - RSHTBR) IS: ",T51,F8. 5) 
57355 FORMAT (' ak SINCE THE NO. OF MISSION TESTS IS THE SAME FOR' 
eter COMPONENTS THE CLOSED FORM SERIES SYSTEM RELIABILITY ' 
‘RHTSTR' IS COMPUTED 
6666 FORMAT (' +! See eee ae ae 
+ | gekdedinticteicine RUN INPUT SETTINGS Jnnididicidinicicicicincicineictenieieinicieteicteteicte 
+ LEAR ERS ER SEEN IC ) 
6667 FORMAT (' ',//'NUMBER OF REPLICATIONS: ',T50,14) 
6668 FORMAT (' ',//'NUMBER OF COMPONENTS: ,T50, 14) 
6669 FORMAT (' ' ue SYSTEM RELIABILITY FUNCTION: ,T50, ‘SERIES' ) 


4447 FORMAT 
4448 FORMAT 
4449 FORMAT 
5555 FORMAT 
5556 FORNAT 


t 
S: 
1 
ae 
S:',T50,F8. 4) 
’ 
4 
1 
' 
? 
t 
5557 FORMAT (‘ ' 
t 
t 


s 
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6699 FORMAT (' ',//'SYSTEM RELIABILITY FUNCTION: ',T50, BRIDGE’) 

6670 FORMAT (' ',//'MASTER UNRELIABILITY USED: ',T50,F8. 5) 

6671 FORMAT (' ',//'INPUT WEIGHTS(A SUB I''S):') 

6674 FORMAT ('+',/// ‘totnieniniiclatrodniiniinntibi old ieedtaivedta det iaksra action rirr | 
+'**R UN RES U_L T Ser ddinioonohntmnininiinininininininticntiricta ! 
4 driededrdictetede voices |) . 

6675 FORMAT (1+! ,/// Gattegroni tee a in rar re nner TaD ISTE TAIN | 
+' ESTIMATE ERRORS * et niidintdddordinintin anda ninbidninntinitdcicicnsnicict | 
+! tokictiniividicict ! 

6676 FORMAT ('+! /// Oittiieaintriinin niet didoinctoc nn inicnicinrn i aincicn icici | 
+' TRUE CONFIDENCE LIMITS eetrrninincineircrniiotdniisciaiinininaiioiiciciverint | 
+ fried decreed dercserea |) 


END 
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APPENDIX B. INTERVAL ESTIMATING PROCEDURE 2 
PROGRAM ZFNSCY . 
kkekkthkhikkkikKitchckicicickktticiccicticnictcktiickiciccdcocekictctcicicsetecickiciicetciciciecletctcicte 
* ¥ 
* TITLE: BINOMIAL INTERVAL ESTIMATION PROCEDURE ve 
* ZERO FAILURES DISALLOWED; WITH SCALING * 
* AUTHOR: E. F. BELLINI, LT, USN * 
* DATE: NOV 89 * 
i ¥e 
* THIS PROGRAM COMPUTES THE TRUE CONFIDENCE LEVEL FOR THE ESTIMATE * 
% RELIABILITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY * 
% OF THEIR COMPONENTS ¥ 
te * 
v IN TTS PRESENT CONFIGURATION THIS PROGRAM IS SET UP TO RUN 12 ve bs 
~o¥ TIMES EACH TIME PRODUCING 1000 REPLICATIONS USING A DIFFERENT ve 
¥ SET OF INPUT DATA. RUN THE PROGRAM FROM CMS BY TYPING ‘Bl EXEC'.* 
ve THE REXX EXEC PROGRAM ¥ 
w "B1i' CALLS THE INPUT FILES TO BE READ AND NAMES THE 12 OUTPUT * 
* FILES RESULTING FROM THE 12 CONSECUTIVE RUNS. BY EDITING THE ¥ 
¥ INDEX COUNTERS I, J, K OF THE 'Bi' EXEC ONE CAN RUN ANY USER- * 
* SPECIFIC RUN FROM JUST ONE RUN TO ALL 12. # 
ve + 
¥ VARIABLES USED * 
ve % 
¥ AHATI : WEIGHT ESTIMATES FOR EACH COMPONENT * 
¥ Al : INPUT WEIGHTS FOR EACH COMPONENT * 
¥ ALFA  : LEVELS OF SIGNIFICANCE * 
* BIGF : TOTAL NO. OF FAILURES FOR EACH REPLICATION * 
* CHISQ : CHI-SQUARE RANDOM VARIABLE VALUE * 
ws COUNT1 : COUNTS THE NO. OF COMPONENTS WITH FAILURES ¥ 
cD C1015 : FORMAT LABEL * 
* DEGFR : DEGREES OF FREEDOM * 
ve DELBRG : DIFFERENCE FOR BRIDGE SYSTEM * 
te DELSTR : DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM * 
ve DELTAR : DIFFERENCE FOR SERIES SYSTEM * 
v DIFF : DIFFERENCE (TRUE REL. - ESTIMATED REL.) | te 
* EPS : SMALL QUANTITY( CONSTANT) * 
% ERROR : PARAMETER FOR IMSL ROUTINE * 
* FAILS : COUNTS NO. OF REPLICATIONS WITH AT LST. 1 FAILURE * 
% FI : NO. OF FAILURES FOR EACH COMPONENT(ALL MISSION TST)* 
xe FLAG : 1 IF ALL COMP. HAVE SAME NO. OF MISSION TESTS * 
* INC : INCREMENT STEP SIZE FOR ROUTINE USMNMX ¥ 
ca KEY1 : ARRAY OF INDECES FOR ROUTINE SHSORT * 
* KEY2  : ARRAY OF INDECES FOR ROUTINE SHSORT * 
* KEY3  : ARRAY OF INDECES FOR ROUTINE SHSORT * 
co KEY4 : ARRAY OF INDECES FOR ROUTINE SHSORT * 
* KK : ARRAY SIZING PARAMETER FOR THE MAX NO OF COMPONENTS* 
¥ LOOP : COUNTS NO. OF REPLICATION PERFORMED 
* MAXALF : MAX NO. OF SIGNIFICANT LEVELS DESIRED( ARRAY sizINey# 
v MAXREP : MAX NO. OF REPLICATIONS 
ve MAXRUN : MAX NO. OF PROGRAM ITERATIONS ALLOWED : 
ve MSTRQ : MASTER UNRELIABILITY(USED WITH Ai'S TO CALC. QI'S) * 
4} 








7 MULT : MULTIPLIER FOR RANDOM NO. GENERATOR SRND 


+ 


* N : NO. OF MISSION TEST FOR EACH COMPONENT * 
* NIMAX : MAX NO. OF MISSION TESTS * 
* NIMIM : MIN NO. OF MISSION TESTS * 
¥e NINDX : INDEX NO. OF MAX NO. OF MISSION TESTS * 
¥ NTREAL : NO. OF MISSION TESTS TRANSFORMED TO REAL * 
* NMAX =: MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL * 
* NPRNT : FLAG FOR DETAILED REPORT OUTPUT * 
* PRNT ; SAME AS ABOVE( PARAMETER) * é 
* QHATI : UNRELIABILITY ESTIMATES FOR EACH COMPONENT * 
* QHTMAX : LARGEST QHATI * 
* QHTUPR : UPPER LIMIT ON SYSTEM UNRELIABILITY * 
* QI : INPUT UNRELIABILIY FOR EACH COMPONENT * 
* QINDX : INDEX * 
* QUANTL : QUANTILE * 
w REPSHD : REPLICATIONS HEADING FORMAT NUMBER * 
* RHTSTR : SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) * 
* RS : TRUE SERIES SYSTEM RELIABILITY * 
* RSBRDG : TRUE BRIDGE SYSTEM RELIABILITY * 
ve RSHAT : SERIES SYSTEM RELIABILITY ESTIMATE * 
* RSHTBR : BRIDGE SYSTEM RELIABILITY ESTIMATE * 
ve SEED : PARAMETER * 
Ye SELCTA : SIGNIFICANCE LEVEL SELECTION * 
* SELCTB : SIGNIFICANCE LEVEL SELECTiC!" ¥ 
* SORT : PARAMETER FOR ROUTINE SRND * 
| ve SUMNAI : SUM OF THE PRODUCT OF NI'S AND AI'S * 
¥ TEMP : TEMPORARY ARRAY * 
| od TOTREP : TOTAL NUMBER OF PROGRAM ITERATIONS * 
* TRANBR : TEMPORARY ARKAY % 
¥ TRANSQ : TEMPORARY ARRAY * 
ve TRANSR : TEMPORARY ARRAY * 
* TRIALS : BERNOULLI TRIALS ARRAY (2-DIM) * 
ve TRNSTR : TEMPORARY ARRAY * 
* TRUQNT : TRUE QUANTILE * : 
* UNIRV +: UNIFORM RANDOM DEVIATES (2-DIN) ve 
w éFAILS : TOTAL NUMBER OF REPLICATIONS WITH ZERO FAILURES * 


ne ZrFPREP : NO. OF COMPNTS. WITH ZERO FAILURES PER REPLICATION * 
* 


Vevetere terete Tee Fete rere vere He We se Ve Me Fe Fe Ve PEL He Fe Le NE LS HTC TE TC He He LET HEHE HELE TE FONE TS He HEHE FEMI Lo HEH LE POE HEI TE NE BE TE ICTS TENE HEI HE 


PARAMETER (KK=10,MAXALF=2 ,NPRNT=0) 
PARAMETER (MAXREP=1000, NAXRUN=10000, EPS=. 000001) 


REAL*4 UNIRV(15, 1000) ,TEMP(1000),QI(KK), AICKK), AHATI(KK) 
REAL*4 QHATI(KK), NMAX, QHTMAX, ALFA(MAXALF) 

REAL*4 DF(5),AALFA(5) ,SUMNAI , RSHAT(MAXALF , MAXREP) ,RS 

REAL*4 KEY1(MAXREP) , KEY2(MAXREP) , KEY3(MAXREP) , TRNSTR(MAXREP) 
REAL*4 DEGFR(MAXREP), QHTUPR(MAXALF ,MAXREP) , CHISQ( MAXALF , MAXREP) 
REAL*4 RHTSTR(MAXALF , MAXREP) 

REAL*4 DELTAR(MAXALF), TRANSQ(MAXREP) , TRANSR(MAXREP) , DIFF(MAXREP) 
REAL*4 DELSTR(MAXALF), NN(KK) 

REAL*4 RSHTBR(MAXALF ,MAXREP) , DELBRG(MAXALF) , KEY4(MAXREP) 

REAL#4 TRANBR(MAXREP), RSBRDG 

REAL*4 SCALEN, MINGI, MAXQI, ZFPREP,MSTRQ 


INTEGER SEED, MULT, SORT, TRIALS(15,1000), BIGF, FI(KK), N(KK) 
















































































INTEGER NINDX, QINDX, ERROR, REPS, SELCTA, SELCTQ, TOTREP 
INTEGER COUNTI, c1c15, REPSHD, SELCTB, ALF, FLAG, LOOP, PRNT 

.. INTEGER QUANTL(MAXALF) , TRUQNT(MAXALF) , 2FAILS, FAILS 
INTEGER INC 


- DATA SEED/123457/, MULT/1/, INC/1/ 
DATA AALFA/.01,.05,.9,.95,.99/, DF/1,5,10,30,40/ 
DATA ALFA/. 20,.050/ 


; * DATA MSTRQ/0.1, .05/ 
* DATA AI/.20,.10,.05,.05,.05/ 
¥ DATA AI/.20,.10,.05,.05,.05,.05,.05,.05,.05,.05/ 
% DATA AI/.05,.05,.05,.05,.05,.05,.05,.05,.05,.05/ 
* DATA AI/.05,.05,.05,.05,.05,.05,.05,.05,.05,.05,.05,.05,.05, 
% -+,05,.05/ 
* DATA-AL/s Upc h Tos V5 Wy 05,405 ,225595,.5/ 


* DATA AI/.1,.1,1.,.05,.2/ 
ig DATA _N/100, 100, 100, 100, 100/,, 


* DATA"N/10,6,5,2,10,4,5,6,7,2/ 

ve DATA N/50,50,50,50,50,50,50,59,50,50/ 

% DATA N/100,100, 100,100,100, 100,100, 100, 100, 100/ 

* DATA N/300,300,300,300,300,300,300, 300,300, 300/ 

% DATA N/20,20,20,20,20,20,20,20,20,20/ 

% DATA N/1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 


i +1000,1000,1000,1000/ 
DATA SORT/0/ 


| ASSIGN 8 TO C1C15 
ASSIGN 9 TO REPSHD 


id CALL EXCMS('FILEDEF 01 DISK B OUTPUT Al (LRECL 132 om 
. % CALL EXCMS('FILEDEF 02 DISK JNK OUTPUT st (LRECL 132 ') 
SELCTQ = 
. PRNT = NPRNT 
DO 12 I=1,KK 
AI(I) = 9999, 


NCI) = 99999999 
12 CONTINUE 
READ( 03 ,*)K,MSTRQ 


DO 11 1=1,K 
READ(03,%*) AICI) ,N(I) 
_ 11. CONTINUE 


vex 1 / INITIALIZE THE QHTUPR ARRAY OF UNRELIABILITY REPLICATIONS, //*#* 





% RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS * 
¥ AND RHTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN * 
week 1 / ALL THE COMPONENT MISSION TESTS ARE EQUAL IN NUMBER [| [ire 


DO 172 ALF=1,MAXALF 
DO 173 REPS=1,MAXREP 
QHTUPR( ALF ,REPS) = 
E RSHAT(ALF ,REPS) = 
RUTSTR(ALF,REPS) = 


id vine sae : Hae 
, sie 




















RSHTBR(ALF,REPS) = 0. 
173 ~ CONTINUE = a “an 
172 CONTINUE 


wie // SET FLAG TO 1 IF ALL COMPONENTS HAVE SAME NO. OF HISSION TESTS*ri#x 


FLAG=1 
DO 50 I=1,K -1 
IFCCN(I) - NCI+1)).NE.0) THEN 
FLAG=0 


ELSE 
END IF 
50 CONTINUE 
PRINT *, 'FLAG IS:', FLAG 


#7 / MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION) //*** 


ZFPREP = 0. 
FAILS = 0 
ZFAILS = 0 
COUNT1 = 0 
TOTREP = 0 
LOOP = 0 


10 IF( LOOP. LT. MAXREP) THEN 
IF(TOTREP. LT. MAXRUN) THEN 
SELCTA = 1 
“<SELCTB = 2 


weve? FILL ARRAY KEY(REPS) WITH INTEGERS 1 TO K TO BE USED AS OUTPUT 


wevere?/ OF THE SUBROUTINE SHSORT 


DO 95 REPS=1, MAXREP 
KEY1(REPS) = REPS 


KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEY4(REPS) = REPS 


95 CONTINUE 


ee othe 


were // CALCULATE NMAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE //*** 


al 


were // CALCULATE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NO. 
CALL IMAX(N,K,NMAX ,NINDX) 


week // CALCULATE THE QI'S FROM THE GIVEN MASTER Q AND THE AI'S 


DO 115 I=1, K 
QI(I) = MSTRQ * AICI) 
115 CONTINUE 


DO 120 I=1,15 
DO 125 J=1,500 
UNIRV(I,J) = 999. 
TRIALS(I,J) = 99999 
125 CONTINUE 
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tae 





ona” 


120 CONTINUE 


wint// CALCULATE THE SCALING FACTOR 'SCALEN' | [itt 


CALL USMNMX(QI,K, INC ,MINQI ,MAXQI) 
SCALEN = MAXQI / MINQI 


wie // DRAW UNIFORM (0,1) RV'S AND CONVERT TO BERNOULLI TRIALS | [*%% 


DO 130 I=1, K 
CALL SRND(SEED, TEMP, N(I), MULT, SORT) 
DO 135 J=1, N(I) 
UNIRV(I,J) = TEMP(J) 
IF (UNIRV(I,J).LE. 1 - QI(I)) THEN 
TRIALS(I,J) = 0 
ELSE 
TRIALS(I,J) 
END IF 


1) 
rary 


135 CONTINUE 
130 CONTINUE 


iek// CALCULATE THE NO. OF FAILURES FOR EACH COMPONENT | [wee 


DO 150 I=1, K 
FI(I) = 0 
150 CONTINUE 


weieve// CALCULATE THE F SUB I'S AND THE GRAND TOTAL NO. OF FAILURES 


BIGF = 0 
DO 155 I=1, K 
DO 160 J=1, N(I) 
FI(I) = FI(I) + TRIALS(I,J) 
160 — m CONTINUE a iad 


a 


wer // CALCULATE THE QHAT SUB I'S: F SUB I'S DIVIDED BY N SUB I'S = //**¥ 
IF(FICI).EQ.0) THEN 
QHATI(I) = 1. / (SCALEN * N(I)) 
ZFPREP = ZFPREP + 1 
SE 
QHATI(I) = REAL(FI(1I)) / NCI) 
IF 


BIGF = BIGF + FI(I) 
155 CONTINUE 


wie // FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES [ [eee 
CALL RMAX(QHATI, K, QHTMAX, QINDX) 


IF(PRNT. EQ. 1) THEN 
WRITE(1,0007) 
WRITE(1,C1C15) 
WRITE(1,3334) QI 
WRITE(1,0001) 











140 


145 


doled / / 


162 





WRITE(1,C1C15) 
DO 140 J=1,NMAX 
WRITE(1,1111) (UNIRV(I,J), I=1, K) 

CONTINUE 

WRITE(1,0002) 

WRITE(1,C1C015) 

DO 145 J=1,NMAX 

WRITE(1,2222) (TRIALS(I,J), I=1, K) 

CONTINUE 

WRITE(1,0003) 

WRITE(1,C1C15) 

WRITE(1,3333) FI 

WRITE(1,0005) 

WRITE(1,C1015) 

WRITE(1,3335) N 

WRITE(1,0004) 

WRITE(1,C1C15) 
...WRITE(1,3334) QHATI 

(1,'C/''THE MAXIMUM Q HAT SUB I IS:'', T40, F8.5)') QHTMAX 
WRITE(1,'(/''THE MAXI Q HAT SUB I IS ELMNT NO.: ' ,T40,15)') QINDX 
WRITE(1,'(/''THE GRAND TOTAL NO. OF FAILURES IS:'',T40, 15)') BIGF 
ELSE 
ENDIF 


TEST FOR A REP WITH AT LST TWO COMP WITH AT LST GNE FAILURE EACH 


DO 162 J=1, K 

NN(J) = 0 

IF(FI(J).GE.1) THEN 
COUNT1 = COUNT1 + 1 
NN(J)= NCJ) ~ 


ELSE | sil 
NN(J) = N(S) * SCALEN 
END IF 
CONTINUE 
IF(COUNT1. EQ. 0) THEN 
ZFAILS = ZFAILS + 1 
ELSE 
FAILS = FAILS + 1 
END IF 
IF(COUNT1.GE. 2) THEN 
LOOP = LOOP + 1 
TOTREP = TOTREP + 1 
ELSE 
TOTREP = TOTREP + 1 
GO TO 10 
END IF 
CALCULATE THE AHAT SUB I'S (WEIGHT ESTIMATES) | [ek 
SUMNAI = 0. 


DO 165 I=1, K 7 
AHATI(I) = QHATI(I) / QHTMAX 
+(FI(L) ,L=1,K) 
+(NN(L) ,L=1,K) 
 SUMNAIL = SUMNAI + NN(I) * AHATI(I) 
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165 CONTINUE 
IF(PRNT. EQ. 1) THEN 
WRITE(1,0006) 
WRITE(1,€1C15) 
WRITE(1,3334) AHATI 
ELSE 


END IF 


wr // CALCULATE 1 REPLICATION OF UPPR ALFA C.L. ON SYSTEM RELIABILITY 
DEGFR(LOOP) = 2 * (1 + BIGF) 


DO 170 ALF=1, MAXALF 
_ CALL MDCHI(1 - ALFACALF) ,DEGFR( LOOP) ,CHISQ(ALF,LOOP), ERROR) 
“QHTUPR(ALF, LOOP) = CHISQ(ALF,LOOP) / (2 * SUMNAT) 
IF(FLAG. EQ. 1) THEN 
eee RHTSTR(ALF,LOOP) = 1 -(CHISQCALF,LOOP) / REAL(2*N(1))) 
END IF 


* + (ALF ,LOOP), ALFACALF) 
wie // CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR COMPNTS. IN SERIES 


CALL RHTSRS(QHTUPR(ALF, LOOP), AHATI,K, RSHAT(ALF, LOOP) ) 
* +T40,F8.5)') RSHAT(ALF, LOOP) 


vei 77 CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE *** 


CALL RHTBRG( QHTUPR( ALF, LOOP) , AHATI ,K,RSHTBR( ALF , LOOP) ) ae! 
170 CONTINUE | 


weet // THIS ELSE AND ENDIF ARE FOR THE TEST AGAINST MAXRUN / [eee 
ELSE 
WRITE(1,'('' '! ,/') PROGRAM EXCEEDED THE MAX NO. OF RUNS'', 
+'' ALLOWED OF: '',16)') TOTREP 
GOTO 9999 


END IF 
GOTO 10 
END IF 
wee 7 / SORT THE ARRAYS OF SYSTEM UNRELIABILITIES(1 FOR EACH CONF. LEVEL) 


DO 700 ALF=1, MAXALF 
DO 800 REPS=1, MAXREP 


TRANSQ(REPS) = QHTUPR(ALF,REPS) 

TRANSR(CREPS) = RSHAT( ALF ,REPS) 

TRNSTRCREPS) = RHTSTRCALF, REPS) 

TRANBRCREPS) = RSHTBRCALF, REPS) 
800 CONTINUE 


CALL SHSORT(TRANSQ, KEY1,MAXREP) 
CALL SHSORT(TRANSR, KEY2 ,MAXREP) 
CALL SHSORT(TRNSTR,KEY3 ,MAXREP ) 
CALL SHSORT(TRANER , KEY4 ,MAXREP) 
DO 900 REPS=1, MAXREP 


47 





QHTUPR(ALF ,REPS) = TRANSQ(REPS) 
RSHAT(ALF,REPS) = TRANSR(REPS) 
RHTSTR(ALF REPS) = TRNSTR(REPS) 
RSHTBR(ALF REPS) = TRANBR(REPS) 

900 |. CONTINUE. i 

700 CONTINUE bi 


wiek // PRINT OUTPUT REPORT HEADINGS 


WRITE( 1, 6666) 
WRITE(1,6667) MAXREP 
WRITE(1,6668) K 
WRITE(1,6770) SCALEN 
WRITE( 1, 6669) 
IF(K.EQ.5) THEN 
WRITE(1,6699) 
ELSE 
END IF 
WRITE(1,6670) MSTRQ 
WRITE(1,6671) 
WRITE(1,C1C15) 
WRITE(1,3334) AI 
WRITE(1,0007) 
WRITE(1,C1C15) 
WRITE(1,3334) QI 
WRITE(1,0005) 
WRITE(1,C1C15) 
WRITE(1,3335) N 
.... WRITE( 1 ,6674) 
WRITE(2, 6666) 
WRITE(2,6667) MAXREP 
WRITE(2,6668) K 
WRITE(2,6770) SCALEN 
WRITE(2,6669) 
IF(K. EQ.5) THEN 
WRITE(1,6699) 
ELSE 
END IF 
WRITE(2,6670) MSTRQ 
WRITE(2,6671) 
WRITE(2,C1C15) 
WRITE(2,3334) AI 2 
WRITE(2,0007) 
WRITE(2,C1C15) 
WRITE(2,3334) QI 
WRITE(2,0005) 
WRITE(2,C1C15) 
WRITE(2,3335) N 
WRITE( 2, 6674) 
WRITE(2,'(''SORTED RSHAT 1 IS:'',/15(F8.5))') 
+(RSHAT(1,REPS), REPS=1, MAXREP) 
WRITE(2,'(''SORTED RSHAT 2 IS:'',/15(F8.5))') 
+(RSHAT(2,REPS), REPS=1, MAXREP) 
IF(FLAG. EQ. 1) THEN 
WRITE(2,'(''SORTED RHTSTR 1 IS:'',/10(F8.5))') 
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+(RHTSTR(1,REPS) , REPS=1, MAXREP) 
WRITE(2,'(° "SORTED RHTSTR 2 IS:'',/10(F8.5))') 
+(RHTSTR(2,REPS), REPS=1, MAXREP) 
ELSE 
END IF 
IF(K. EQ.5) THEN 
WRITE(2,'(''SORTED RSHTBR 1 IS:'',/10(F8.5))') 
+(RSHTBR(1,REPS) , REPS=1, MAXREP) 
WRITE(2,'(° "SORTED RSHTBR 2 IS:'',/10(F8.5))') 
+(RSHTBR(2,REPS), REPS=1, MAXREP) 
ELSE 
END IF = na ase 





wee // AND FOR THE 5-COMPONENT BRIDGE STRUCTURE 


CALL RSRS(QI,K,RS) 
WRITE(1,'('' '',////''THE TRUE SERIES SYSTEM '', 
+'' RELIABILITY VALUE IS: '',T51,F8.5)') RS 

CALL RBRIDG(QI ,K,RSBRDG) 

; IFCK. EQ. 5) THEN 


WRITE(1,’('' '',////''THE TRUE BRIDGE STRUCTURE '', 
; +'' RELIABILITY VALUE 1S:'',T51,F8.5)') RSBRDG 

ELSE 

END IF 


WRITE(1,6675) 


weere// COMPUTE THE DIFFERENCE 'DELTAR' BTWN. RS AND RSHAT OF THE THEO 
wet // RETICAL QUANTILE GIVEN BY ALFA(MUST USE SORTED RSHAT ARRAY) 


IF(FLAG. EQ. 1) THEN 
WRITE(1,5755) 
LSE 
5 END IF 
DO 450 ALF=1, MAXALF 
QUANTL(ALF) = MAXREP * (1 - ALFA(ALF)) 
DELTAR(ALF) = RS - RSHAT(ALF,QUANTL(ALF)) 
DELBRG(ALF) = RSBRDG - RSHTBR( ALF, QUANTL( ALF) ) 
IF(FLAG. EQ. 1) THEN 
DELSTR(ALF) = RS - RHTSTR(ALF,QUANTL(ALF)) . 
WRITE(1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 
! WRITE(1,5656) RHTSTR(ALF , QUANTL( ALF) ) 
WRITE(1,5657) DELSTR(ALF) 
ELSE © 
END IF 
IF(K.EQ.5) THEN 
DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTL(ALF) ) 
WRITE(1,5555) MAXREP, ALFAC(ALF), REAL(QUANTL(ALF)) 
WRITE(1,5666) RSHTBR(ALF ,QUANTL( ALF) ) 
WRITE(1,5667) DELBRG( ALF) 
ELSE 
END IF 
™  WRITE(1,5555) MAXREP, ALFA(ALF) ** REAL(QUANTL( ALF) ) 
. WRITE(1,5556) RSHAT(ALF , QUANTL(ALF) ) 
. WRITE(1,5557) DELTAR(ALF) 
450 CONTINUE 





*k%// COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL. FNCTN. (SERIES SYSTEM) 














PRINT *, ‘ QUANTL( 1) IS:', QUANTL(1) ERIS ae. 


PRINT x, *QUANTL(2) IS:', QUANTL(2) 


#77 FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*¥** 


wick / / tidied: RSHAT teeeineser / [ee 
WRITE(1,6676) ~ “ ™ ee 7 ee 
DO 400 ALF=1,MAXALF 
TRUQNT(ALF) = 


DO 500 REPS=1, MAXREP 
DIFF(REPS) = RS - RSHAT(ALF,REPS) 
500 CONTINUE 
DO 600 REPS=1, MAXREP 
IF(ABS(DIFF(REPS)). ne tee THEN 


TRUQNT(ALF) = 
WRITE(1,'('' ue oy TRUE CONFIDENCE LIMIT Is:'' 
+ F8.4)') 
+ CTRUQNTC ALF) / REAL(MAXREP)) * 100. 
GO TO 620 


ELSEIF(DIFF(REPS). LT. 0. ) THEN 
TRUQNT( ALF) = REPS 


GO TO 610 
ELSE 
END IF 
600 CONTINUE 
610 IF( TRUQNT( ALF). EQ. 0.) THEN 


WRITE(1,4443) ALFA( ALF) 
WRITE(1,'('' '',/''THE SMALLEST'', 
+ "' DIFFERENCE BETWEEN RS AND RSHAT IS:'',F10.5)') DIFF( 
+ MAXREP) 
ELSEIF(TRUQNT( ALF). EQ. 1.) THEN 
WRITE(1,4442) ALFA( ALF) 
WRITE(1,'('' '',/'' ALL RSHAT'’ 


op "' ARE GREATER THAN RS'’)' Na? 
ELSEIF( ABS( DIFF(TRUQNT( ALF) )). LE. ABS(CDIFF(TRUQNT(ALF) - 1))) 
+ THEN 
WRITE(1,4444) ALFA(ALF), 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 


WRITE(1,4445) RSHAT( ALF, TRUQNT( ALF) ) 
WRITE(1,4446) 
ELSE 
WRITE(1,4444) ALFA(ALF), 
+ ((TRUQNT(ALF)-1) / REALC(MAXREP)) * 100. 
WRITE(1,4445) RSHAT( ALF , TRUQNT(ALF)-1) 
WRITE(1,4447) 
620 END IF 
400 CONTINUE 


sine » ea! cote 





week // FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*** 
aevede J f wetetedeseveteve RSHTBR (BRIDGE) *ieteeetve [lick 


IF(K. EQ. 5) THEN 
DO 401 ALF=1,MAXALF 
TRUQNT( ALF) = 0 

DO 501 REPS=1, MAXREP 


50 











DIFFC(REPS) = RSBRDG - RSHTBR(ALF,REPS) 
501 CONTINUE 


DO 601 REPS=1, MAXREP 
IF(ABS(DIFF(REPS)). LE. EPS) THEN 
p TRUQNT(ALF) = REPS 
: WRITE(1, "('' '',/''TRUE CONFIDENCE LIMIT IS:'', 


+ F8.4)') 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
GO TO 621 


ELSEIF(DIFF(REPS). LT. 0.) THEN 
TRUQNT( ALF) = REPS 


GO TO 611 
ELSE 
END IF 
601 CONTINUE 
611 IF(TRUQNT( ALF). EQ. 0.) THEN 


WRITE(1,4443) ALFA( ALF) 
WRITE(1,' ('' ''\7''THE SMALLEST'', 
eile My DIFFERENCE BETWEEN RSBRDG AND RSHTBR IS: Le os 
+ F10.5)') DIFF(MAXREP) 
ELSEIF(TRUQNT( ALF). EQ. 1.) THEN 
WRITE(1, 4442) ALFA( ALF) 
WRITE(1,' Sa eeeanre T RSHTBR'', 

+ ‘' ARE GREATER THAN RSBRDG'')' ) 

ELSEIF( ABS( DIFF(TRUQNT( ALF) )). LE. ABS(DIFF(TRUQNT(ALF) - 1))) 

+ THEN 
WRITE(1,4444) ALFA(ALF), 

+ (TRUQNT( ALF) / REAL(MAXREP)) * 100. 
WRITE(1,4449) RSHTBR( ALF , TRUQNT( ALF) ) 
WRITE(1,4446) 

ELSE 
coe WRITE(1,4444) ALFACALF), wy 
+ _((TRUQNT(ALF)-1) / REAL(MAXREP)) * 100. 
. " WRITE(1,4449) RSHTBR( ALF , TRUQNT( ALF) -1) 
WRITE(1,4447) 
621 END IF 
401 CONTINUE 
ELSE 
END IF 


wiek 77 FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //%** 
tededs / / Kicdedstetetele RHUTSTR wetervededevededs / / teiede 


IF(FLAG. EQ. 1) THEN 
DO 4400 ALF=1,MAXALF 
TRUQNT(ALF) = 0 
DO 5500 REPS=1, MAXREP 
DIFF(REPS) = RS - RHTSTR(ALF,REPS) 
5500 CONTINUE 
DO 6600 REPS=1, MAXREP 
IF(ABS(DIFF(REPS)). LE. EPS) THEN 
TRUQNT(ALF) = REPS 
WRITE(1,'('' '',/''TRUE CONFIDENCE LIMIT IS: '' 
+ Fs. 4)') ' 
+ (TRUQGNT(ALF) / REAL(MAXREP)) * 100. 
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60 TO 6620 ai 
ELSEIF(DIFF(REPS). LT. 0.) THEN 
TRUQNT( ALF) = REPS 


GO TO 6610 
ELSE 
END IF 
6600 CONTINUE 
6610 IF(TRUQNT( ALF). EQ. 0.) THEN 


WRITE(1,4443) ALFA(ALF) 
WRITE(1,'('' °',/''THE SMALLEST'', 
+ { DIFFERENCE BETWEEN RS AND RHTSTR IS:'' 
+ F9.5)') DIFF(MAXREP) 
ELSEIF(TRUQNT( ALF). EQ. 1.) THEN 
WRITE(1,4442) ALFA( ALF) 
WRITE(1,' ('' "',/''ALL RHTSTR'', 


+ "" ARE GREATER THAN RS'')"') 
a ee ae net LE. ABS(DIFF(TRUQNT(ALF) - 1))) 
+ HE 
WRITE(1,4444) ALFACALF), 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
WRITE(1,4448) RHTSTR(ALF , TRUQNT( ALF) ) 
WRITE(1,4446) © 
ELSE 
WRITE(1,4444) ALFACALF), 
+ (C(TRUQNT(ALF)-1) / REAL(MAXREP)) * 100. 


WRITE(1,4448) RHTSTRCALF , TRUQNT( ALF) -1) 
wp WRITE(1,4447) 
6620 END IF 
4400 CONTINUE 
ELSE 
END IF 


yee / 7 PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION //*#* 


IF(PRNT. EQ. 1) THEN 
185 WRITE(1,REPSHD) ALFACSELCTA), ALFA(SELCTA), 
+ALFA(SELCTB) ,ALFAC(SELCTB) , ALFA(SELCTA) , ALFACSELCTA) , ALFAC(SELCTB), 
+ALFACSELCTB) 
175. IFCI.GE. (MAXREP + 1)) THEN 
“GOTO 180 a “ ne 
ELSE He 
IF( (I. EQ. 71). OR. (I. EQ. 211). OR. (I. EQ. 351). OR. (I. EQ. 491). OR. 
+ (1. EQ. 631). OR. (1. EQ. 771). OR. (I. EQ. 911). OR. (I. EQ. 1051) ) THEN 
I= I+ 70 ; 
WRITE(1,'(''+'")") 
GOTO 185 
ELSE 


WRITE(1,3336) I, INT(DEGFR(I)), CHISQ(1,1), QHTUPR(1,I), 
+ CHISQ(2,1), QHTUPR(2,1) 
END IF 
IF(I + 70. LE. MAXREP) THEN 
WRITE(1,3337) 1+70, INT(DEGFR(1+70)) ,CHISQ(1,I+70), 
+ QHTUPR(1,I+70) ,CHISQ(2,1+70) ,QHTUPR(2,1+70) 
ELSE . 


Pry 


LD 
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END IF ae 
I=I+1 
GOTO 175 
180 END IF 
ELSE 
ENDIF 
9999 WRITE(1,'(''THE TOTAL NO OF REPS WAS:'',I8)') TOTREP 
WRITE(1,'(''THE TOTAL NO OF EFFECTIVE REPS WAS:'' 18)") LOOP 
WRITE(1,'(''THE TOTAL NO OF NO FAILURE RUNS was: '! 18) ) ZFAILS 
WRITE(1,'('' THE TOTAL NO OF RUNS WITH FAILURES WAS: '',1I8)') FAILS 
Ra ae OF COMPONENTS WITH NO FAILURES PER REP WAS:'', 
. a {AXRE 
0008 FORMAT (/ 3X,'C 1' ,5X,'C of, 
#5X%-'C 3) 5X, C6" {5X,°6. 5 5X, C6" 5X5 °C. 7" 5k. 
+'C 8" 5X,'C 9',5X,'C 10',4X,'C 11',4x, 
+'C 12° ,4X,'C 13° ,4X,'C 14" ,4X,'C 15') 
0009 FORMAT(/1X, 'REP NO" ,2X,'DF »1X,'CHISQR(' ,F4.3,')",1X 
+'QHTUPR(' ,F4.3,')',1X, CHISQR( ;F4.3,')" 1X, 'QHTUPR( »F4.3,')', 
+2X,'REP NO’ ,2X,'DF » 1X," CHISQR( »F4.3,' ,1X, 
+'QHTUPR(' ,F4.3,')',1X, ‘CHISQR(' ,F4. 3, )' ,1X, QHTUPR(',F4.3,')'/) 
0001 FORMAT (///'UNIFORM RANDOM DEVIATES ARE: ) 
0002 FORMAT (///'BERNOULLI TRIALS ARE: ') - 


0003, FORMAT (///'TOTAL NO. OF FAILURES FOR EACH COMPONENT: ') 
~ 0004 FORMAT (///'ESTIMATED UNRELIABILITY FOR EACH COMPONENT: ') 


0005 FORMAT 


(///' TOTAL NUMBER OF MISSION TESTS: ') 

0006 FORMAT (/ 
(/ 

(1 


//'ESTIMATED WEIGHTS FOR EACH COMPONENT: ') 
//'Q I FOR EACH COMPONENT: ') 
5F8. 5) 


0007 FORMAT 
1111 FORMAT 


2222 FORMAT (/1X,15(14,4X)) 
3333 FORMAT (/1X,15(14,4X)) 
3334 FORMAT (/15F8. 5) 
3335 FORMAT (/1X,15(14,4X)) 
3336 FORMAT (T3,14,T9,13,T13,F11.5,T27,F8.5,T39,F11.5,T53,F8. 5) 
3337 FORMAT ('+',T67,14,T73,13,T77,F11.5,T91,F8.5,7103,F11.5,T117,F8. 5) 
4442 FORMAT (' ',///'THE RESULTING (1 - ',F4.3,') CONFIDENCE ', 
+'LIMIT IS:',T50,' 00.000. ‘) ae sa 
4443 FORMAT (' ',///'THE RESULTING (1 - ',F4.3,') CONFIDENCE ', 
+' LIMIT ",T50,'100. 0000') 


( 
( 
( 
( 
( 
( 
( 
I 
( 
IS: 
4444 FORMAT (' ' ///' THE RESULTING (1 - ',F4.3,') CONFIDENCE ', 
+'LIMIT IS: ',T50,F8. 4) 
4445 FORMAT (' 
4446 FORMAT (' 
4447 FORMAT (' 
4448 FORMAT (' 
4449 FORMAT (' 
5555 FORMAT (' 
5556 FORMAT (' 
5557 FORMAT (' 
5656 FORMAT (' 
5666 FORMAT (' 
5657 FORMAT (' 
5667 FORMAT (' 


',/'THE RSHAT VALUE CLOSEST TO RS IS: ',T51,F8.5) 
',/'CFIRST NEGATIVE DIFFERENCE)" ) 

',/'CELEMENT PRECEEDING FIRST NEGATIVE DIFFERENCE)‘ ) 
',/'THE RHTSTR VALUE CLOSEST TO RS IS: ',T51,F8.5) 
',/'THE RSHTBR VALUE CLOSEST TO RSBRDG IS: ',T51,F8.5) 
",///'THE ',14,'(1-',F4.3,') QUANTILE IS: ',T49,F8. 3) 
',/'THE VALUE OF RSHAT FOR THAT QUANTILE IS: ',T51,F8.5) 
',/'THE DIFFERENCE(RS - RSHAT) IS:',T51,F8.5) 

',/'THE VALUE OF RHTSTR FOR THAT QUANTILE IS: ',T51,F8.5) 
',/'THE VALUE OF RSHTBR FOR THAT QUANTILE IS:',T51,F8.5) 
'\/'THE DIFFERENCE(RS - RHTSTR) IS: ',T51,F8.5) 

',/'THE DIFFERENCE(RS - RSHTBR) IS:',T51,F8.5) 


5755 FORMAT (' ',///'SINCE THE NO. OF MISSION TESTS IS THE SAME FOR', 
+' ALL COMPONENTS THE CLOSED FORM SERIES SYSTEM RELIABILITY ', 
+'''RHTSTR'' IS COMPUTED' ) 


in 
ne 














6666 FORMAT (1+! Gitininn orn toininiicgantaiinintoi iia igoi a i aii iain a irieicie a | 


+'* RUN IS IPUT SETTINGS FOR ADJUSTED SUM OF THE NI AI  eefeticiininiict | 
+! Pootseabdere de aikrese de akstatserc tes ! ) 


6667 FORMAT (' ',//'NUMBER OF REPLICATIONS: ',T50,14) 


6668 FORMAT (' ',//'NUMBER OF COMPONENTS: ',T50,14) 

6669 FORMAT (' ',//'SYSTEM RELIABILITY FUNCTION: ',T50,'SERIES') 
6699 FORMAT (' ',//'SYSTEM RELIABILITY FUNCTION: ',T50, ' BRIDGE’ ) 
6670 FORMAT (' ',//'MASTER UNRELIABILITY USED: ,T50,F8. 5) 

6770 FORMAT es ZERO FAILURE SCALING FACTOR: ° »T50,F5. 2) 


6671 FORMAT ,//'INPUT WEIGHTS(A SUB I''S): 

6674 FORMAT ('+' ,///'‘iotdddnainidicoinigoidioigicoincinaiciotddi iat aigiigicotitctoiiat: | 
+'%*R UN RE S U LT Sttnininteiohitninidninineeiidiiioiddiciietiticta | 
+ ieteiiiniinticitn |) 

6675 FORMAT (14! ,/// Uitisddinineiicct i tiation ata ic Ir I TRI ICI RT TOR RITA ARI 
+' ESTIMATE ERRORS (}tikiniigdinieioiinnidctnniciieininicnidiinincicicinincn inti: | 
4+ doldintoticintiatictiy | ) 

6676 FORMAT ('+!' | /// Unnigniddinididiniioiioidiodnidninigininicietnncotictdincinincicicticcn | 
+' TRUE CONFIDENCE LIMITS tte triniedeinidinetedidiotinigcaicniotincicisnticia i! 
+! HMekokvek eek arar a | ) 


END 


$4 








APPENDIX C. INTERVAL ESTIMATING PROCEDURE 3 


PROGRAM ZFYSCY 
Hebert Perec AR ITI TRE A ITT RETR IIR TAIT NR TCR TTR I IIIT TA II TRI TA NAAT A IN 


* * 
* TITLE: BINOMIAL INTERVAL ESTIMATION PROCEDURE * 
* ZERO FAILURES DISALLOWED; WITH SCALING * 
* AUTHOR: E. F. BELLINI, LT, USN * 
* DATE: NOV 89 * 
* * 
* THIS PROGRAM COMPUTES THE TRUE CONFIDENCE LEVEL FOR THE ESTIMATE * 
* RELIABILITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY * 
* OF THEIR COMPONENTS ve 
+ * 
ve IN ITS PRESENT CONFIGURATION THIS PROGRAM IS SET UP TO RUN 12 * 
* TIMES EACH TIME PRODUCING 1000 REPLICATIONS USING A DIFFERENT ¥ 
* SET OF INPUT DATA. RUN THE PROGRAM FROM CNS BY TYPING ‘bl EXEC'.* 
¥ THE REXX EXEC PROGRAM ve 
* "B1' CALLS THE INPUT FILES TO BE READ AND NAMES THE 12 OUTPUT ve 
7 FILES RESULTING FROM THE 12 CONSECUTIVE RUNS. BY EDITING ThE ¥ 
* INDEX COUNTERS I, J, K OF THL 'B1' EXEC ONE CAN RUN ANY USER- ¥ 
SPECIFIC RUN FROM JUST ONE RUN TO ALL 12. vr 
5 id bid 
% VARIABLES USED 7 
big * 
¥ AHATI : WEIGHT ESTIMATES FOR EACH COMPONENT w 
ve Al : INPUT WEIGHTS FOR LACH COMPONENT * 
° a ALFA : LEVELS OF SIGNIFICANCE * 
oo BIGF : TOTAL NO. OF FAILURES FOR EACH REPLICATION a 
¥ CHISQ : CHI-SQUARE RANDOM VARIABLE VALUE * 
% COUNT1 : COUNTS THE NO. OF COMPONENTS WITH FAILURES * 
* C1C15  : FORMAT LABEI. * 
bag DEGFR : DEGREES OF FREEDOM 7 
7 DELBRKG : DIFFERENCE FOR BRIDGE SYSTEM sd 
DELSTR : DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM * 

* LELTAR : DIFFERENCE FOR SERIES SYSTEM * 
bs DIFF : DIFFERENCE (TRUE REL. - ESTIMATED REL. ) ve 
i EPS : SMALL QUANTITY ( CONSTANT ) ¥ 
ERROR : PARAMETER FOR IMSL ROUTINE 7 

* FAILS : COUNTS NO. ot} REPLICATIONS WITH Al LST. 1 FAILURE * 
* FI : NO. OF FAILURES FOR EACH COMPONENTCALL MISSION TST)* 
= FLAG > 1 JF ALL Cui’. HAVE SAME NO. OF MISSION TESTS * 
% INC : INCREMENT STEP SIZE FOR ROUTINE USMNMX - 
* KEY1 : ARRAY OF INDECES FOR ROUTINE SHSORT * 
a KEY2 : ARRAY OF INDECES FOR ROUTINE SHSOKT * 
* KLY3 : ARRAY OF INDECES FOR ROUTINE SHSORT * 
as KEYS : ARRAY OF INDECES FOR ROUTINE SHSORT ve 
* KK : ARRAY SIZING PARAMETER FOR THE MAX NO OF COMPONENTS* 
* LOOP : COUNTS NO. CF REPLICATIGN PERFORMED ne 
% MASALF : SAX NO. OF SIGNIFICANT LEVELS DESIREDC ARRAY SiZING)* 
7 MASREP : MAX NO. OF REPLICATIONS * 
* * MAXRUN : MAX NO. OF PROGRAM ITCRATIONS ALLOWED * 


a MSTRQ : MASTER UNKELIABLLITY(USED WITH AL'S TO CaLC. QI'S) * 





+ 


MULT : MULTIPLIER FOR RANDOM NO. GENERATOR SRND 
N : NO. OF MISSION TEST FOR EACH COMPONENT * 
NIMAX : MAX NO. OF MISSION TESTS 
NIMIM : MIN NO. OF MISSION TESTS 
NINDX : INDEX NO. OF MAX NO. OF MISSION TESTS 
NIREAL : NO. OF MISSION TESTS TRANSFORMED TO REAL 
NMAX =: MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL 
NPRNT : FLAG FOR DETAILED REPORT OUTPUT 
PRNT  : SAME AS ABOVE( PARAMETER) 
: UNRELIABILITY ESTIMATES FOR EACH COMPONENT 
QHTMAX : LARGEST QHATI 
QHTUPR : UPPER LIMIT ON SYSTEM UNRELIABILITY 
Qi : INPUT UNRELIABILIY FOR EACH COMPONENT 
QINDX : INDEX 
QUANTL : QUANTILE 
REPSHD : REPLICATIONS HEADING FORMAT NUMBEk 
RHTSTR : SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) 
RS : TRUE SERIES SYSTEM RELIABILITY 
RSBRDG : TRUE BRIDGE SYSTEM RELIABILITY 
RSHAT : SERIES SYSTEM RELIABILITY ESTIMATE 
RSHTBR : BRIDGE SYSTEM RELIABILITY ESTIMATE 
SEED : PARAMETER 


+etetee ee etter t t+ $44 4+ + 
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SELCTA : SIGNIFICANCE LEVEL SELECTION * 
* SELCTB : SIGNIFICANCE LEVEL SELECTION * 
* SORT : PARAMETER FOR ROUTINE SRND * 
* SUMNAI : SUM OF THE PRODUCT OF NI'S AND AI'S * 
i TEMP : TEMPORARY ARRAY * 
. TOTREP : TOTAL NUMBER OF PROGRAM ITERATIONS * 
* TRANBR : TEMPORARY ARRAY * 
* TRANSQ : TEMPORARY ARRAY * s 
* TRANSRK : TEMPORARY ARRAY ad 
* TRIALS : BERNOULLI TRIALS ARRAY (2-DIM) = 
* TRNSTR : TEMPORARY ARRAY * 
ve TRUQNT : TRUE QUANTILE * 
* UNIRV : UNIFORM RANDOM DEVIATES (2-DIM) * 
* ZFAILS : TOTAL NUMBER OF REPLICATIONS WITH ZERO FAILURES * 
1 ZFPREP : NO. OF COMPNTS. WITH ZERO FAILURES PER REPLICATION * 
* * 


Hlededevedede te detetesevedeseh sete ve Tete sete te dese te te deve de te te eye Fe Fe Ve Ve Ve Ne He Te Fe Fe Fe LF HC TC He Te He FE HE LC Le TETHER LEKI RI HE 


PARAMETER (KK=10,MAXALF=2 ,NPRNT=0 ) 
PARAMETER (MAXREP=1000, MAXRUN=10000, EPS=. 000001) 


REAL*4 UNIRV( 15,1000), TEMP( 1000) ,QI(KK), AI(KK), AHATICKK) 
REAL*4 QHATI( KK), NMAX, NNMAX, QHTMAX, ALFA(MAXALF ) 
REAL*4 DF(5),AALFA(5) ,SUMNAI ,RSHAT(MAXALF ,MAXREP) ,RS 
REAL*4 KEY1(MAXREP) ,KEY2(MAXREP) , KEY3( MAXREP) , TRNSTR( MAXREP ) 
REAL*4 DEGFR(MAXREP), QHTUPR( MAXALF ,MAXREP) , CHISQ(MAXALF , MAXREP ) 
REAL*4 RHTSTR( MAXALF ,MAXREP ) 
REAL*4 DELTAR(MAXALF), TRANSQ(MAXREP) , TRANSR(MAXREP) , DIFF(MAXRE? ) 
REAL*4 DELSTR(MAXALF), NN( KK) 
REAL*4 RSHTBR(MAXALF ,MAXREP) , DELBRG(MAXALF) , KEYS( MAXREP ) 
REAL*4 TRANBRC(MNAXREP), RSBRDG 
REAL*4 SCALEN, MINQI, MAXQI, ZFPREP,MSTRQ 


INTEGER SEED, MULT, SORT, TRIALS(15,10C0), BIGF, FICKK), NCKK) 









































INTEGER NINDX, QINDX, ERROR, REPS, SELCTA, SELCTQ, TOTREP 
INTEGER COUNTI, c1C15, REPSHD, SELCTB, ALF, FLAG, LOOP, PRNT 
INTEGER QUANTL(MAXALF) , TRUQNT(MAXALF) , ZFAILS, FAILS 
INTEGER INC 


DATA SEED/123457/, MULT/1/, INC/1/ 
DATA AALFA/.01,.05,.9,.95,.99/, DF/1,5,10,30,40/ 
DATA ALFA/. 20,.050/ 


* DATA MSTRQ/0.1, .05/ 

Xt DATA AI/.20,.10,.05,.05,.05/ 

x DATA AI/.20,.10,.05,.05,.05,.05,.05,.05,.05,.05/ 
% DATA AI/.05,.05,.05,.05,.05,.05,.05,.05,.05,.05/ 


% DATA AI/.05,.05,.05,.05,.05,.05,.05,.05,.05,.05,.05,.05,.05, 
% +,05,.05/ 
% DATA AI/.1,.1,.1,.1,1,-05,.05,.2,.25,.5/ 


% DATA AI/.1,.1,1.,.05,.2/ 

* DATA N/100,100, 100,100, 100/ 

* DATA N/10,6,5,2,10,4,5,6,7,2/ 

* DATA N/50,50,50,50,50,50,50,50,50,50/ 

x DATA N/100,100,100,100,100,100, 100,100,100, 100/ 

** DATA N/300, 300,300, 300,300,300, 300, 300,300, 300/ 

* DATA N/20,20,20,20,20,20,20,20,20,20/ 

* DATA N/1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000 , 1000, 1000, 





+1000,1000,1000, 1000/ 
DATA SORT/0/ 


ASSIGN 8 TO C1C15 
ASSIGN 9 TO REPSHD 





. ve CALL EXCMS('FILEDEF 01 DISK B OUTPUT Al (LRECL 132 Ds 

a CALL EXCMS(' FILEDEF 02 DISK JNK OUTPUT J Al (LRECL 132 ? 
SELCTQ = 
PRNT = NPRNT 


DO 12 I=1,KK 
AICI) = 9999. 
NCI) = 99999999 
12 CONTINUE 
READ( 03 ,*)K,MSTRQ 


DO 11 [=1,K 
READ(03,*) AICI) ,N(1) 
11 CONTINUE 


wick // INITIALIZE THE QHTUPR ARRAY OF UNRELIABILITY REPLICATIONS, [ [eee 
* 


RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS * 
% AND RHTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN * 
vekte// ALL THE COMPONENT MISSION TESTS ARE EQUAL IN NUMBER [ [2% 


DO 172 ALF=1,MAXALF 
DO 173 REPS=1,MAXREP 
QHTUPR(ALF,REPS) = 
. RSHAT(ALF,REPS) = 
RHTSTR(ALF,REPS) = 0. 


57 








ree SS 


RSHTBR(ALF,REPS) = 0. 
173 CONTINUE 
172 CONTINUE 


wiek// SET FLAG TO 1 IF ALL COMPONENTS HAVE SAME NO. OF MISSION TESTS#ie% 


FLAG=1 
DO 50 I=1,K -1 
IF((N(I) - NCI+1)).NE.0) THEN 
FLAG=0 


ELSE 
END IF 
50 CONTINUE 
_..PRINT *, ‘FLAG IS:', FLAG 
wie 77° MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION) //*** 


ZFPREP = 

FAILS = 0 

ZFAILS 

COUNT1 

TOTREP 

LOOP = 0 

10 IF(LOOP. LT. MAXREP) THEN 

IF(TOTREP. LT. MAXRUN) THEN 

SELCTA = 1 

SELCTB = 2 


0 
0 
0 


0 


wrk // FILL ARRAY KEY(REPS) WITH INTEGERS 1 0 X TO BE USED .AS OUTPUT 
#i%// OF THE SUBROUTINE SHSORT 


DO 95 REPS=1, MAXREP 


KEY1(REPS) = REPS 
KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEY4(REPS) = REPS 


95 CONTINUE | | 

wee // CALCULATE NMAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE 

rr// CALCULATE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NO. | / [ee 
CALL IMAX(N,K,NMAX,NINDX) 


yt] / CALCULATE THE QI'S FROM THE GIVEN MASTER Q AND THE AI'S | [tee 


DO 115 I=1, K 
QI(I) = MSTRQ * AI(I) 
115 CONTINUE 





DO 120 I=1,15 
DO 125 J=1,500 
UNIRV(I,J) = 999. 
TRIALS(1,J) = 99999 
125 CONTINUE 
120 CONTINUE 





any 
w 


wint// CALCULATE THE SCALING FACTOR 'SCALEN’ [ [wie 


CALL USMNMX(QI,K, INC ,MINQI ,MAXQI) 
SCALEN = MAXQI / MINQI 


DRAW UNIFORM (0,1) RV'S AND CONVERT TO BERNOULLI TRIALS //** 


. DO 130 I=1, K 
CALL SRND(SEED, TEMP, N(I), MULT, SORT) 
DO 135 J=1, NCI) 
UNIRV(T,J) = TEMP(J) 
IF CUNIRV(I,J).LE. 1 - QI(1I)) THEN 
0 


: Yetere / / 


TRIALS(I,J) = 
ELSE 
TRIALS(I,J) = 1 
END IF 
135 CONTINUE 
130 CONTINUE 
wk%// CALCULATE THE NQ. OF FAILURES FOR EACH COMPONENT (ES = 


DO 150 I=1, K 
FI(I) = 
150 CONTINUE 


wiv / / CALCULATE THE F SUB I'S AND THE GRAND TOTAL NO. OF FAILURES 


BIGF = 0 
- DO 155 I=1, K 
DO 160 J=1, N(I) 
FICI) = FICI) + TRIALS(I,J) 
160 CONTINUE 


wee / 7 CALCULATE THE QHAT SUB I'S: F SUB I'S DIVIDED BY N SUB I'S 
IF(FI(1I).£Q.0) THEN 
QHATI(I) = 1. / (SCALEN * N(I)) 
ZFPREP = ZFPREP + 1 
ELSE . 
QHATI(I) = REAL(FI(1I)) / NCI) 
END IF 
BIGF = BIGF + FICI) 
155 CONTINUE 





ert FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES 
CALL RMAXCOHATI, K, QHTMAX, QINDX) 


IF(PRNT. EQ. 1) THEN 
WRITE(1,0007) 
. WRITE(1,C1C15) 
, WRITE(1,3334) QI 
WRITE(1,0001) 
WRITE(1,C1C15) 
DO 140 J=1,NMAX 
- WRITE(1,1111) (UNIRV(I,J), I=1, K) 
140 CONTINUE | 








~. WRITE(1,0002) 
WRITE(1,C1C15) 
DO 145 J=1,NMAX 
WRITE(1,2222) (TRIALS(I,J), I=1, K) 
145 CONTINUE 
WRITE(1,0003) 
WRITE(1,C1C15) 
WRITE(1,3333) FI 
WRITE(1,0005) 
WRITE(1,C1C15) 
WRITE(1,3335) N 
WRITE(1,0004) 
WRITE(1,C1C15) 
~MRITE( 1, 3334) QHATI ~ 
WRITE(1,"(/' THE MAXIMUM Q HAT SUB I Ts: » THO, FB. 5)') QHTMAX 
WRITE(1,'(/''THE MAXI Q HAT SUB I IS ELMNT NO.: ',T40,15)") QUNDX 
WRITE(1,'(/''THE GRAND TOTAL NO. OF FAILURES Is:'!,T40, 15)! ) BIGF 
ELSE 
ENDIF 


wit // TEST FOR A REP WITH AT LST TWO COMP WITH AT LST ONE FAILURE EACH 
DO 162 J=1, K 
NN(J) = 0 
IF(FI(J).GE.1) THEN 
COUNT1 = COUNT1 + 1 


NN(J)= NCJ) 
ELSE 

NN(J) = NCJ) * SCALEN 
END IF 


162 CONTINUE ° 
IF(COUNT1. EQ. 0) THEN 
ZFAILS = ZFAILS + 1 
ELSE 
FAILS = FAILS + 1 
END IF 
IF(COUNT1.GE.0) THEN 
LOOP = LOOP + 1 
TOTREP = TOTREP + 1 
ELSE 
TOTREP = TOTREP + 1 
GO TO 10 
END IF 





week 1/7 CALCULATE THE AHAT SUB I'S (WEIGHT ESTIMATES) 


SUMNAI = 0. 
DO 165 I=1, K 
AHATICI) = QHATI(I) / QHTMAX 
SUMNAI = SUMNAI + NN(I) * AHATI(I) 
165 CONTINUE 
IF(PRNT. EQ. 1) THEN 
~ WRITE( 1, 0006) ae me _ 
WRITE(1,C1C15) . 
WRITE(1,3334) AHATI 
ELSE 
END IF 
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wie // CALCULATE 1 REPLICATION OF UPPR ALFA C.L. ON SYSTEM RELIABILITY 
DEGFR(LOOP) = 2 * (1 + BIGF) 


DO 170 ALF=1, MAXALF 


CALL MDCHI(1 - ALFA(ALF) ,DEGFR(LOOP) ,CHISQ(ALF,LOOP), ERROR) 
QHTUPR(ALF,LOOP) = CHISQ(ALF,LOOP) / (2 * SUMNAT) 
_IF(FLAG. EQ. 1) THEN 

: RHTSTRCALF , LOOPY = 1 -(CHISQCALF,LOOP) / REAL(2*N(1))) .°™ 
ELS 
END IF 


wiek// CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR COMPNTS. IN SERIES. - 
CALL RHTSRS(QHTUPR(ALF,LOOP), AHATI,K, RSHAT(ALF, LOOP) ) 
* +T40,F8.5)') RSHAT(ALF, LOOP) 
week // CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE *** 
CALL RHTBRG( QHTUPR( ALF , LOOP) , AHATI ,K,RSHTBR( ALF , LOOP) ) 
170 CONTINUE 
yee // THIS ELSE AND ENDIF ARE FOR THE TEST AGAINST MAXRUN [ [we 
ELSE 
WRITE(1,'(' '! /'! PROGRAM EXCEEDED THE MAX NO. OF RUNS'', 
+'' ALLOWED OF: '' a6 ) TOTREP 
GOTO 9999 
END IF 
GOTO 10 
END IF 
wee SORT THE ARRAYS OF SYSTEM UNRELIABILITIES(1 FOR EACH CONF. LEVEL) 


DO 700 ALF=1, MAXALF 
DO 800 REPS=1, MAXREP 


TRANSQ(REPS) = QHTUPR(ALF ,REPS) 

TRANSR(REPS) = RSHAT(ALF,REPS) 

TRNSTR(REPS) = RHTSTR(ALF,REPS) 

TRANBR(REPS) = RSHTBR(CALF,REPS) 
600 CONTINUE 


CALL SHSORT(TRANSQ, KEY1,MAXREP) 

CALL SHSORT( TRANSR , KEY2 ,MAXREP) 

CALL SHSORT(TRNSTR, KEY3 ,MAXREP) 

_CALL SHSORT( TRANBR ,KEY4,MAXREP) — 

“DO 900 REPS=1, MAXREP 8s a 
QHTUPR(ALF ,REPS) = TRANSQ(REPS) 
RSHAT( ALF ,REPS) = TRANSR(REPS) 
RHTSTR(ALF ,REPS) = TRNSTR(REPS) 
RSHTBR(ALF ,REPS) = TRANBR( REPS) 

990 CONTINUE 
700 CONTINUE 


6] 








wor // PRINT OUTPUT REPORT HEADINGS / [ick 


WRITE(1,6666) 
WRITE(1,6667) MAXREP 
WRITE(1,6668) K« sk 
WRITE(1,6770) SCALEN a , 
WRITE(1,6669) 

IF(K.EQ.5) THEN 
WRITE(1, 6699) 
ELSE 
END IF 
WRITE(1,6670) MSTRQ 
WRITE(1,6671) 
WRITE(1,C1C15) 
WRITE(1,3334) AI 
WRITE(1,0007) 
WRITE(1,C1C15) 
WRITE(1,3334) QI 
WRITE(1, 0005) 
WRITE(1,01C15) 
WRITE(1,3335) N 
WRITE(1,6674) 
WRITE(2,6666) 
WRITE(2,6667) MAXREP 
WRITE(2,6668) K 
WRITE(2,6770) SCALEN 
WRITE(2, 6669) 
IF(K.EQ.5) THEN 
WRITE(1, 6699) 

.. ELSE 

END IF’ 

WRITE(2,6670) MSTRQ 

WRITE(2,6671) 

WRITE(2,C1C15) 

WRITE(2,3334) AI 

WRITE( 2 ,0007) 

WRITE(2,C1C15) 

WRITE(2,3334) QI 

WRITE(2,0005) 

WRITE(2,C1C15) 

WRITE(2,3335) N 

WRITE( 2 ,6674) 

WRITE(2,'(''SORTED RSHAT.1 IS:'',/15(F8.5))') 
+(RSHAT(1,REPS), REPS=1, MAXREP) 
WRITE(2,'(''SORTED RSHAT 2 IS:'',/15(F8.5))") 
+(RSHAT(2,REPS), REPS=1, MAXREP) 

IF(FLAG. EQ. 1) THEN 

WRITE(2,'(''SORTED RHTSTR 1 IS:'',/10(F8.5))') 
+(RHTSTR(1,REPS), REPS=1, MAXREP) 
WRITE(2,'(''SORTED RHTSTR 2 IS:'',/10(F8.5))') 
+(RHTSTR(2,REPS), REPS=1, MAXREP) 
ELSE 
END IF 
IF(K.EQ.5) THEN 
WRITE(2,'(''SORTED RSHTBR 1 IS:'',/10(F8.5))') 
+(RSHTBR(1,REPS), REPS=1, MAXREP) 





+ PIR 
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~WRITE(2,'(“SORTED RSHTBR 2 IS:+',/10(F8.5))') 
+(RSHTBR(2,REPS), REPS=1, MAXREP) 
ELSE 
END IF 


aay 


wrk // COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL. FNCTN. (SERIES SYSTEM) 
week // AND FOR THE 5-COMPONENT BRIDGE STRUCTURE 


ad CALL RSRS(QI,K,RS) 
WRITE(1,'C('' '',////''THE TRUE SERIES SYSTEM ea 
+' "RELIABILITY VALUE IS:'' »T51,F8.5)') RS 
CALL RBRIDG(QI,K,RSBRDG) 
IF(K. EQ. 5) THEN 


WRITE(1,'('' '',////''THE TRUE BRIDGE STRUCTURE ue 
+'' RELIABILITY VALUE Is: '' »T51,F8.5)') RSBRDG - 

ELSE 

END IF 


WRITE(1,6675) 


wee / 7 COMPUTE THE DIFFERENCE 'DELTAR' BTWN. RS AND RSHAT OF THE THEO 
wie // RETICAL QUANTILE GIVEN BY ALFA(MUST USE SORTED RSHAT ARRAY) 


IF( FLAG. EQ. 1) THEN 
WRITE(1,5755) 


ELSE 
’ END IF 
DO 450 ALF=1, MAXALF 
QUANTL( ALF) = MAXREP * (1 - ALFACALF)) 
DELTARCALF) = RS - RSHATCALF ,QUANTL( ALF) ) 
. DELBRG(ALF) = RSBRDG - RSHTBR( ALF, QUANTL( ALF) ) 


IF(FLAG. EQ. 1) THEN 
DELSTR(ALF) = RS - RHTSTR(ALF,QUANTL(ALF)) 
, WRITE(1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 
WRITE(1,5656) RHTSTR(ALF, QUANTL( ALF) ) 
WRITE(1,5657) DELSTR(ALF) 
ELSE 
END IF 
IF(K.EQ.5) THEN 
vn DELBR6( ALF) =-RSBRDG - RSHTBR(ALF,QUANTL(ALF) ) a 
- WRITE(1,5555) MAXREP, ALFA(ALF), REAL( QUANTL( ALF) ) | 
WRITE(1,5666) RSHTBR( ALF ,QUANTL( ALF) ) 
WRITE(1,5667) DELBRG(ALF) 
ELSE 
END IF 
WRITE(1,5555) MAXREP, ALFACALF), REAL(QUANTL(ALF)) 
WRITE(1,5556) RSHAT(ALF, QUANTL(ALF)) 
WRITE(1,5557) DELTAR(ALF) 
450 CONTINUE 


*x*// FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*** 
= kek / / Fertekiciciekke RSHAT veteacdedeaede / / teicte 
WRITE( 1, 6676) 

DO 400 ALF=1,MAXALF 
TRUQNT( ALF) = 0 
DO 500 REPS=1, MAXREP 
DIFFCREPS) = RS - RSHAT(ALF,REPS) 
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500 CONTINUE 
DO 600 REPS=1, MAXREP 
IF(ABS(DIFF(REPS)).LE.EPS) THEN 
TRUQNT( ALF) = REPS 


WRITE(1,'('' '',/''TRUE CONFIDENCE LIMIT IS:'', 
+ F8.4)') i 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

GO TO 620 


ELSEIF(DIFF(REPS). LT. 0.) THEN 
TRUQNT(ALF) = REPS 


GO TO 610 
ELSE 
END IF 
600 CONTINUE 
610 IF(TRUQNT( ALF). EQ: 0.) THEN 
WRITE(1, 4443) ALFA( ALF) 
WRITE(1,'('' '',/''THE SMALLEST'’, 
+ '' DIFFERENCE BETWEEN RS AND RSHAT Is: '" »F10.5)') DIFF( 
+ MAXREP) 


ELSEIF(TRUQNT(ALF). EQ. 1.) THEN 
WRITE(1,4442) ALFA(ALF) 


: WRITE(1,'('' '',/''ALL RSHAT'', 
+ '' ARE GREATER THAN RS'')") 
ELSEIF(ABS(DIFF(TRUQNT( ALF) ).). LE. ABS(DIFF(TRUQNT(ALF) - 1))) 
+ THEN 
WRITE(1,4444) ALFA(ALF), 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 


WRITE(1,4445) RSHAT(ALF , TRUQNT( ALF) ) 
WRITE(1,4446) 
ELSE 
WRITE(1,4444) ALFA(ALF), 
+ ((TRUQNT(ALF)-1) / REAL(MAXREP)) * 100. 
WRITE(1,4445) RSHAT( ALF, TRUQNT(ALF)-1) 
WRITE(1,4447) 
620 END IF ~ sie om 
&00 CONTINUE 


“7 / FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE | [ik 
wet set adoekekiek RSHTBR £ (BRIDGE) detevedededede [ [ere 


IF(K. EQ. 5) THEN 
DO 401 ALF=1,MAXALF 
TRUQNT(ALF) = 0 
DO 501 REPS=1, MAXREP 
DIFF(REPS) = RSBRDG - RSHTBR(ALF,REPS) 
501 CONTINUE 
DO 61 REPS=1, MAXREP 
IFC ABS(DIFF(REPS)).LE.EPS) THEN 
TRUQNT( ALF) = REPS 


WRITE(1,'("' '' | /''TRUE CONFIDENCE LIMIT IS:'', 
+ F8.%)') 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

GO TO 621 


ELSEIF(DIFFCREPS). LT. 0.) THEN 
TRUQNT( ALF) = REPS 
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GO TO 611 
ELSE 
END IF 
601 CONTINUE 
611 IF(TRUQNT( ALF). EQ. 0.) THEN 


WRITE(1,4443) ALFA(ALF) 

WRITE(1,'('' '',/"'THE SMALLEST'', 
+ '' DIFFERENCE BETWEEN RSBRDG AND RSHTBR IS:'', 
+ F10.5)') DIFF(MAXREP) 

ELSEIF(TRUQNT( ALF). EQ. 1.) THEN 

WRITE(1,4442) ALFA( ALF) 

WRITE(1,'('' '',/''ALL RSHTBR'', 
+ '' ARE GREATER THAN RSBRDG''’)') 
ELSEIF( ABS(DIFF(TRUQNT( ALF) )). LE. ABS(DIFF(TRUQNT(ALF) - 1))) 
THEN 


WRITE(1,4444) ALFACALF), 
a Eb Pan (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
WRITE(1,4449) RSHTBR( ALF, TRUQNT( ALF) ) 
WRITE(1,4446) 
ELSE 
WRITE(1,4444) ALFA(ALF), 
+ ((TRUQNT(ALF)-1) / REAL(MAXREP)) * 100. 
WRITE(1,4449) RSHTBR( ALF, TRUQNT( ALF) -1) 
WRITE(1,4447) 
621 END IF 
401 CONTINUE 
ELSE 
END IF 


vment oe ate pepe 


sete] 1 FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*** ~ 
vereve / rf Fekietertevetek RHTSTR wredededcdeaede / / Fetete 


cect tet 


IF(FLAG. EQ. 1) THEN 
DO 4400 ALF=1,MAXALF 
TRUQNT(ALF) = 0 
DO 5500 REPS=1, MAXREP 
DIFF(REPS) = RS - RHTSTR(ALF,REPS) 
5500 CONTINUE 
DO 6600 REPS=1, MAXREP 
IF(ABS(DIFF(REPS) ). LE. EPS) THEN 
TRUQNT(ALF) = REPS 
WRITE(1,'('' '',/''TRUE CONFIDENCE LIMIT IS:'', 


+ F8.4)') 
+ (TRUQNT( ALF) / REAL(MAXREP)) * 100. 
GO TO 6620 


ELSEIF(DIFF(REPS). LT. 0. ) THEN 
TRUQNT(ALF) = REPS 
GO TO 6610 
ELSE 
END IF 
€600 CONTINUE 
6610 IF(TRUQNT( ALF). EQ. 0.) THEN 
WRITE(1,4443) ALFA(ALF) 
WRITE(1,'C'' '',/''THE SMALLEST'', - 
+ '' DIFFERENCE BETWEEN RS AND RHTSTR IS:'', . 
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F9.5)') DIFF(MAXREP) 
ELSEIF(TRUQNT( ALF). EQ. 1.) THEN 

WRITE(1,4442) ALFA( ALF) 

WRITEC TC"? tJ" MALL RHTSTR'', 


+ "ARE GREATER THAN RS'')') 
: ELSEIF(ABS( DIFF (TRUQNT( ALF) ) ). LE. ABS(DIFF(TRUQNT(ALF) - 1))) 
WRITE(1,4444) ALFA(ALF), 
+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 
WRITE(1,4448) RHTSTR(ALF, TRUQNT( ALF) ) 
WRITE(1,4446) 
ELSE 
WRITE(1,4444) ALFA(ALF), 
+ ((TRUQNT(ALF)-1) / REAL(MAXREP)) * 100. 
WRITE(1,4448) RHTSTR(ALF, TRUQNTCALF) 1) 
WRITE(1,4447) 
6620 END IF 
4400 CONTINUE 
ELSE 
END IF 


Weve? 7 PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION / /7* 


IF(PRNT. EQ. 1) THEN 


I=1 


Coad ong 
ied ae) 


185 WRITE(1,REPSHD) ALFA(SELCTA), ALFA(SELCTA), 
+ALFA(SELCTB) , ALFA(SELCTB) , ALFA(SELCTA) , ALFA(SELCTA) , SEPA eerere ys 
+ALFA(SELCTB) 


175 


+ 


$999 


TPCT. 


ELSE 


Ils 


GE. (MAXREP + 1)) THEN 
GOTO 180 


IFC (I.EQ. 71). OR. (1. EQ. 211). OR. (1. EQ. 351). OR. (1. EQ. 491). OR. 
(I. EQ. 631). OR. (1. EQ. 771). OR. (1. EQ. 911). OR. (1. EQ. 1051) }) THEN 
I= 1+ 70 
WRITE(1,'(''+'")') 
GOTO 185 
ELSE 


WRITE(1,3336) I, INT(DEGFR(I)), CHISQ(1,1), QHTUPR(1,I), 

CHISQ(2,1), QHTUPR(2,I) 

END IF 

IF(I + 70. LE.MAXREP) THEN 
WRITE(1,3337) I+70.INT(DEGFR(I+70)), CHISQ(1, I+70), 
QHTUPR( 1, I+70) , CHISQ(2,1+70) , QHTUPR( 2, 1+70) 

ELSE | 

END IF. 

I+1 


GOTO 175 


END 


IF 


ELSE 

ENDIF 

WRITE(1,'(''THE TOTAL NO OF REPS WAS: '',I8)') TOTREP 
WRITE(1,'(''THE TOTAL NO OF EFFECTIVE REPS WAS: '',18)') LOOP 
WRITE(1,'(''THE TOTAL NO OF NO FAILURE RUNS WAS: '',I8)') ZFAILS— 
WRITE(1,'("'THE TOTAL NO OF RUNS WITH FAILURES WAS: '',18)') FAILS 
WRITE(1,"(''AVG NO OF COMPONENTS WITH NO FAILURES PER REP WAS:'', 


oway 











475.2)') 2FPREP / MAXREP 


0C08 FORMAT (/ 3X 


"CA" 5X. 2" 


#5, Co 9" OX, te 4! Re to 5! nie ROU) Piste ae) a 
#°C 8° 5K, C2" .5X,"C 10° 4X, C. 11',4X, 
+'C 12° ,4X, 'C 13" 4X, ‘C 14" 4x, "C 15") 
0009 FORMAT( /1X, 'REP No' 2x, 'DF! 1k, CHISQR(' ,F4.3,")",1X 


+' QHTUPR( ' »F4. 3,! )' ; ; 
+2X, "REP NO! 32X, ‘DF y1X " CHISQR( 


ix, CHISQR( 


Fa, 3 a 1X, *outuBRe! »F4.3,')', 


+'QHTUPR(', FA, 3,')', 1x, 'cHTSQRC' fae 3, 3 ik, TQHTUPR(' ,F4.3,')'/) 
ry 


0001 FORMAT 





C7" UNIFORM RANDOM DEVIATES ine: 


0002 FORMAT (///'BERNOULLI TRIALS ARE: ') 
0003 FORMAT (///'TOTAL NO. OF FAILURES FOR EACH COMPONENT: oe 
0004 FORMAT (///'ESTIMATED UNRELIABILITY FOR EACH COMPONENT: ") 
0005 FORMAT (///'TOTAL NUMBER OF MISSION TESTS: ') 
0006 FORMAT (///'ESTIMATED WEIGHTS FOR EACH COMPONENT: ') 
0007 FORMAT (///'Q I FOR EACH COMPONENT: ') 
1111 FORMAT (15F8.5) 
| 2222 FORMAT (/1X,15(14,4X)) 
3333 FORMAT™(/1X,15(14;4X)) nan vm 
3334 FORMAT (/15F8.5) 
3335 FORMAT (/1X,15(14,4X)) 
3336 FORMAT (T3,14,T9,13,T13,F11.5,T27,F8.5,T39,F11.5,1T53,F8. 5) 
3337 FORMAT (' + 767,14, 173,13,T77,F11. $4,731; F8.5, T103,Fil. 5, T117,F8.5) 
4442 FORMAT (' ',/// THE RESULTING (1 - | F4.3,") CONFIDENCE ', 
+'LIMIT IS:',T50,' 00.000 ') 
4443 FORMAT c TL TE RESULTING (1 - ',F4.3,')} CONFIDENCE ', 
+'LIMIT IS: ',T50,'100. 0000’) 
4444 FORMAT (' tp THE RESULTING (1 - ',F4.3,') CONFIDENCE ' 
+'LIMIT IS:',T50,F8. 4) 
4445 FORMAT (' ',/'THE RSHAT VALUE CLOSEST TO RS IS: ';T51,F8.5) 
4446 FORMAT (' ',/'(FIRST NEGATIVE DIFFERENCE)" ) 
4447 FORMAT (' ',/'(ELEMENT PRECEEDING FIRST NEGATIVE DIFFERENCE)') 
4448 FORMAT (' ',/'THE RHTSTR VALUE CLOSEST TO RS IS: ',T51,F8.5) 
4449 FORMAT (' ',/'THE RSHTBR VALUE CLOSEST TO RSBRDG IS: ',T51,F8.5) 
5555 FORMAT (' ',///'THE ',14,'(i-',F4.3,') QUANTILE IS:',T49,F8. 3) 
5556 FORMAT (' ',/'THE VALUE OF RSHAT FOR THAT QUANTILE IS:',T51,F8.5) 
5557 FORMAT (' ',/'THE DIFFERENCE(RS - RSHAT) IS:',T51,F8.5) 
5656 FORMAT (' ',/'THE VALUE OF RHTSTR FOR THAT QUANTILE IS:',T51,F8.5) 
5666 FORMAT (' ',/'THE VALUE OF RSHTBR FOR THAT QUANTILE IS: ',T51,F8.5) 
5657 FORMAT (' vo /\ THE DIFFERENCE(RS - RHTSTR) IS: ',T51,F8.5) 


5755 FORMAT 


6666 FORMAT 





6667 FORMAT 
6668 FORMAT 
6669 FORMAT 
7 6699 FORMAT 
6670 FORMAT 
6770 FORMAT 
‘é 6671 FORMAT 
6674 FORMAT 


- 


5667 FORMAT ( 


/'THE DIFFERENCE(RS - RSHTBR) IS: ',T51,F8.5) 


C O/T ' SINCE THE NO. OF MISSION TESTS IS THE SAME FOR', 
+' ALL COMPONENTS THE CLOSED FORM SERIES SYSTEM ‘RELIABILITY ' 
+''"RHTSTR'' IS COMPUTED') 


, 


('+', " dededededes deeiciciidiefdniciscinniniciicicideiciiccitededsetiedctieniesetye | ‘ 


. ,//' NUMBER OF REPLICATIONS: ' ,T50,14) 
(' ',//'NUMBER OF COMPONENTS: ',T50,14) 


'\//'ZERO FAILURE SCALING FACTOR: ' ,T50,F5. 2) 


i] 
! 
',//'MASTER UNRELIABILITY USED: ' ,TS0,F8. 5) 
7 »//' INPUT WEIGHTS(A SUB I''S):') 


( 
( 
( 
¢ 
( 
CF 


os om Td 


PO 


',//'SYSTEM RELIABILITY FUNCTION: ' ,T50,'SERIES' 
',//' SYSTEM RELIABILITY FUNCTION: ' ,T50,'BRIDGE' 


+'%* RUN INPUT SETTINGS FOR ADJUSTED SUM OF THE NI ATI: tetseierniediicte! 


4! deiededeteiciieteietenniickieict |) 


) 
) 


101] / | Uieiesesicieieinieieieineinteietainieieininicietseiniteieinieininieietciieioinoleicice 


ia 





Ae ne te 


+'**#R UN RESULT Sr Ter eA eA Re IRIE IR TER TNT RH THIN PRT ITE I 
+ incideddinitcatints ! 

6675 FORMAT ('+!' 7 /// otgcdbchiniidadiqhiniciinidotoigd ol tet tnicnnnnkaicincn reir ! 
+' ESTIMATE ERRORS FAIRER III IIT TI RT IR RCT IIR FOR FORTE 
4+ nedosninticiiticint |) 

6676 FORMAT ('+' ,///' htehigidnieddoiniodinicinindinntodincicinniaincrnnorin st | 
+' TRUE CONFIDENCE LIMITS 3etteterdnccibinotatidhiiidntanet tretinoin , 
+ lobinidideodiesticd | 5 


END ‘ 
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APPENDIX D. EXTERNAL SUBROUTINES 


These six subroutines are used by all three programs listed in appendices one 
through three. They must be appended to the particular main program being run or they 
must be available on the same disk as the one from which the main program is being run. 

SUBROUTINE IMAX(SEQ, N, MX, INDEX) 


al THIS ‘ROUTINE IDENTIFIES THE MAXIMUM ELEMENT OF AN INTEGER VECTOR 
NOTE THAT SINCE THE TEST IS .GT. SUBROUTINE ONLY PICKS THE FIRST 


* OCCURRENCE OF A MAX SUCH AS IN THE CASE OF A TIE. 
werede J | ALSO, ARRAY PASSED MUST BE TYPE INTEGER 
REAL*4 MX 
INTEGER SEQ(N), N, I, INDEX 
MX = 0. 
DO 5 I=1, N 
IF(SEQ(1I).GT.MX) THEN 
MX = SEQ(I) 
INDEX = I 
ELSE 
END IF 
5 CONTINUE fi _ 
END 


SUBROUTINE RBRIDG(QI,N,RRSS) 


were // THIS SUBROUTINE CALCULATES THE "TRUE" RELIABILITY OF A 5-COMPONENT 
wwe // BRIDGE STRUCTURE 
REAL*4 QI(N), RRSS 
INTEGER N 
IF(N.NE.5) THEN 
WRITE(1,'(''WARNING: BRIDGE STRUCTURE ONLY USES '', 
+''THE FIRST 5 COMPONENTS'')') 
ELSE 
END IF 
RRSS=( 1-QI(1))*(1-QI(4) )+(1-QI(2) )*(1-QI(5))+(1- Qt(1))*(1- QI(3))* 
CC 1-QI(S) )+C1-QIC2) )**C1-QIC3) )*C1-QI( 4) ) -(1-QI( 1) )*(1-QI(2))* 
C(1-QI(3) )*(1-QI(4) )-€1-QIC 1) )*C1-QIC2) )*(1-QI(3) )*(1-QI(S))- 
CC 1-QIC 1) )*(1-QI(2) )*C1-QI(4) )*(1-QI(5) )-€1-QI( 1) )*(1-QI(3) )* 
CC 1-QI(4) )*(1-QI(5))-€1-QI(2) )*C1-QI(3) )*C1-QI(4) )*(1-QI(5) + 
Nee ty aan eRe BIE eed 
N = 


SUBROUTINE RHTBRG( QHTUP ,AHT,N,RRBRDG) 
wie // THIS SUBROUTINE CALCULATES THE ESTIMATED RELIABILITY OF A 
wie? 1 S-COMPONENT BRIDGE STRUCTURE. (ONLY CARRIED OUT TO THE Q-CUBED TERM 
REAL*4 QHTUP, RRBRDG, AHT(N) 
INTEGER N 
RRBRDG=1-( (QHTUP**2)*( AHT( 1)*AHT(2)+AHT(4)*AHT(5)))-~ 
CC (QHTUP***3 )**( AHT( 1) *ANT(3)*AHT(5S )+AHT( 2) *AHT(3)*ANTC4) ) )+ 
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CC (QHTUP**4 )*( AHT( 1) *AHT( 2) *AHT(3)*AHT( 4)+AHT(1)*AHT( 2) *AHT(3)* 
CAHT(5)+AHT( 1)*AHT( 2)*AHT( 4 )*AHT( 5 )+AHT( 1)*AHT( 3) *AHT( 4) *AHT( 5) + 
CAHT(2)*AHT( 3)*AHT(4)*AHT(5)))- 

C2*( (QHTUP**5 )*( AHT( 1)*AHT( 2) "AHT( 3) *AHT(4)*ANT(5))) 

END 


SUBROUTINE RHTSRS( QHTUP, AAHTI ,N,RRSHAT) 


whe // THIS ROUTINE CALCULATES THE VALUE OF THE SYSTEM RELIABILITY OF A 
* SERIES SYSTEM OF 'N' NO. OF COMPONENTS WHICH HAVE UNRELIABILITY 
wir // 'QETUP'. THE FINAL SYSTEM RELIABILITY VALUE PASSED IS 'RRSHAT' 

REAL*4 QHTUP, RRSHAT, AAHTI(N) 

INTEGER I, N 

RRSHAT = 1. 

no 100 I=1, N 

RRSHAT = RRSHAT * (1 -(CAAHTI(1)* QHTUP)) 
100 CONTINUE 
END 


SUBROUTINE RMAX(SEQ, N, MX, INDEX) 
veer // THID ROUTINE IDENTIFIES THE MAXIMUM ELEMENT OF A REAL VECTOR 
ha NOTE THAT SINCE THE TEST IS .GST. SUBROUTINE ONLY PICKS THE FIRST 
* OCCURRENCE OF A MAX SUCH AS IN THE CASE OF A TIE. 
woe f / ALSO, ARRAY PASSED MUST BE ‘rYPE REAL 
REAL*4 MX, SEQ(N) 
INTEGER N, I, INDEX 


NX = 0. 
DO 5S I=1, N ’ 
IF(SEQ(1).GT.MX) THEN 
MX = SEQ(I) 
INDEX = I 
ELSE , 
END IF 
5 CONTINUE 
END 


SUBROUTINE RSRS(QIS,N,RRS) 


wiei// THIS ROUTINE CALCULATES THE VALUE OF THE SYSTEM RELIABILITY OF A 
ve SERIES SYSTEM OF 'n' COMPONENTS WHICH HAVE UNRELIABILITY 
wae / 'QIS'. THE FINAL SYSTEM RELIABILITY VALUE PASSED IS 'RRS' 
REAL*4 QIS(N), RRS 
INTEGER I, N 
RRS = 1, 
DO 100 I=1, N 
RRS = RRS * (1 - QIS(I)) 
100 CONTINUE 
END 
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